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ABSTRACT 


This  study  addresses  the  problem  of  finding  an  energy  conversion  method  to  take 
advantage  of  the  high  maximum  cycle  temperatures  achieved  with  solar  central  receivers. 
Most  current  practice  is  to  use  steam-based  heat  engines  with  solar  receivers.  Solar 
central  receiver  applications  offer  maximum  cycle  temperatures  in  excess  of  the  operating 
range  of  steam-based  Rankine  cycles.  Gas-based  Brayton  cycles  are  considered  to  take 
advantage  of  the  expected  increased  cycle  performance  of  higher  temperatures. 

Simple  Brayton  cycles  have  been  eliminated  from  consideration  since  they  do  not 
offer  efficiency  improvements  over  the  lower  temperature  Rankine  cycles.  Derivatives 
of  Brayton  engines  do  offer  favorable  performance  compared  to  Rankine  cycles. 
Computer  modelling  was  done  to  examine  the  effect  of  maximum  temperature  on 
efficiency  of  two  Brayton  cycle  derivatives  and  a  simple  Brayton  cycle  (denoted  as  GT). 
The  modified  Brayton  cycles  include  a  combination  of  regeneration  and  one  stage  each  of 
intercooling  and  reheat  (IGT)  and  a  Brayton  cycle  with  steam  injection  (STIG).  The 
turbine  inlet  temperature,  the  steam-to-air  injection  mass  ratio  (for  the  STIG)  and  the 
compression  pressure  ratios  were  treated  as  parameters  in  the  analysis.  The  compression 
ratios  investigated  are  4,  8,  12,  16,  and  20.  Turbine  inlet  temperatures  ranged  from 
1000  K  (1340'F)  to  2500  K  (4040°F).  For  the  STIG  the  steam/air  mass  ratio  varied 
between  0  and  0.5.  Both  First  Law  and  Second  Law  efficiencies  were  examined  as 
turbine  inlet  temperature  and  steam/air  mass  ratios  were  varied  for  each  compression 
ratio. 

Results  show  that  efficiencies  were  highest  for  the  IGT  followed  by  the  STIG  and 
GT,  respectively.  Considerable  improvements  in  specific  work  output  were 
demonstrated  by  the  STIG  over  both  the  IGT  and  GT  systems.  First  Law  and  Second 


Law  analyses  of  the  effect  of  turbine  inlet  temperature  show  a  gradual  increase  of 
efficiency  with  turbine  inlet  temperature  with  diminishing  returns  at  higher  temperatures. 
These  results  demonstrate  reason  for  optimism  about  increased  performance  for  higher 
temperature  air-based  systems  compared  to  Rankine  systems  at  presently  used 
temperatures. 
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INTRODUCTION 


The  current  practice  for  power  generation  from  solar  central  receivers  is  to  use 
steam-based  heat  engines.  Even  the  newer  receiver  concepts,  using  molten  salt  or  liquid 
sodium,  are  envisioned  to  use  a  heat  exchanger  external  to  the  receiver  to  generate  steam 
for  driving  a  Rankine  cycle  system  to  produce  power.  In  spite  of  this  apparent 
contentment  with  steam-based  systems,  it  is  well-known  that  there  should  be  a 
thermodynamic  basis  for  expecting  increasing  cycle  performance  if  the  cycle  maximum 
temperature  is  increased.  However,  any  substantial  increase  in  cycle  maximum 
temperature  may  well  preclude  the  use  of  steam.  Gas-based  cycles  are  therefore 
considered  for  this  application. 

Air  cycles  have  been  considered  for  solar  applications  over  the  years.  The  Stirling 
cycle  and  Brayton  cycle  are  both  undergoing  some  development  for  solar  power 
generation  systems  based  upon  dish  (distributed)  receivers.  For  larger  systems  such  as 
might  be  used  for  solar  central  receivers,  the  development  of  Brayton-derivative  cycles 
seems  to  be  one  of  the  few  alternatives  available.  Brayton  applications  have  been 
generally  eliminated  from  serious  consideration  because  they  appear  to  offer  too  much  of 
an  efficiency  penalty  compared  to  lower  temperature  systems.  This  impression  has  been 
drawn  because  of  increased  receiver  losses  as  temperatures  are  increased  as  well  as  the 
seemingly  low  cycle  performance  of  simple  Brayton  cycles.  However,  there  are  some 
Brayton  engine  derivatives  that  may  demonstrate  favorable  performance  compared  to 
Rankine  systems.  This  study  will  evaluate  the  performance  of  two  of  these  derivatives 
and  compare  them  to  the  performance  of  the  simple  Brayton  cycle  (denoted  as  GT). 


These  include  a  familiar  combination  of  a  Brayton  cycle  with  regeneration  and  one  stage 
each  of  intercooling  and  reheat  (denoted  as  IGT)  and  a  less  investigated  Brayton  cycle 
with  steam  injection  (STIG).  This  study  determines  realistic  performance  of  these  cycles 
operating  on  air.  Special  concern  is  given  to  the  impact  of  maximum  cycle  temperature 
on  overall  performance. 

As  explained  earlier,  current  solar  power  generation  systems  use  a  heat  exchanger 
to  transfer  the  heat  from  the  molten  salt  or  liquid  sodium  fluid  to  a  steam  based  Rankine 
cycle.  Part  of  the  problem  of  making  full  use  of  higher  temperatures  afforded  by  solar 
systems  is  transfer  of  the  solar  energy  to  the  power  generation  system.  This  study  will 
not  address  the  problem  of  how  the  exchange  of  heat  from  the  initial  fluid  medium  of  the 
solar  control  receiver  is  made  to  the  air  or  steam/air  mixture  working  fluid  of  the  cycles 
evaluated.  This  problem  is  a  significant  one  and  is  not  addressed  here.  In  this  study  this 
heat  exchanger  is  termed  as  a  heater  in  the  text  and  Figure  1.  It  will  be  assumed  that  the 
solar  central  receiver  exists  to  heat  the  working  fluid  to  the  turbine  inlet  temperatures 
investigated  in  this  study. 

The  configurations  of  these  three  systems  are  illustrated  in  Figure  1.  The 
advantages  the  IGT  system  enjoys  compared  to  the  simple  GT  system  include  the 
intercooler  to  decrease  compressor  work  and  the  reheat  to  increase  turbine  work.  Since 
both  of  these  components  incur  an  efficiency  penalty,  the  IGT  system  includes  a 
regenerator  to  increase  efficiency  by  preheating  the  compressor  discharge  with  the 
turbine  exhaust.  The  performance  characteristics  of  this  Brayton  derivative  cycle  have 
been  well-documented.  To  serve  as  an  example,  several  gas  turbine  texts  were  reviewed 
to  compare  the  improvements  offered  by  the  IGT  cycle  over  the  GT  cycle.  Boyce  [1] 
reports  an  81.4%  increase  in  specific  net  work  and  a  53.5%  increase  in  efficiency  at  a 
turbine  inlet  temperature  of  1256  K  (1800°F)  and  a  pressure  ratio  of  10.5  At  a  higher 
temperature  of  1478  K  (2200'F)  and  the  same  pressure  ratio,  he  reports  an  increase  of 
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52.5%  in  specific  net  work  and  a  56.2%  increase  in  efficiency.  Bathie  [2]  states  that  the 
IGT  cycle  provides  a  33.3%  increase  in  specific  net  work  and  a  42.2%  increase  in 
efficiency  over  the  GT  cycle  at  a  turbine  inlet  temperature  of  1400  K  (2060T).  Wilson 
[3]  in  his  recent  text  lists  performance  curves  for  a  Brayton  cycle  with  intercooling  and 
regeneration  (no  reheat).  This  system  offered  an  increase  of  20.1%  in  specific  net  work 
and  an  increase  of  23.1%  in  efficiency  over  the  GT  cycle  at  a  turbine  inlet  temperature  of 
1465  K  (2177°F)  and  a  pressure  ratio  of  20.  At  a  higher  turbine  inlet  temperature  of 
2051  K  (3232°F)  and  the  same  pressure  ratio,  he  reports  an  increase  of  30.3%  in 
specific  net  work  and  41.4%  in  efficiency  over  the  GT  cycle. 

The  STIG  system  also  increases  cycle  efficiency  by  use  of  the  turbine  exhaust  to 
generate  steam.  This  steam  is  then  injected  into  the  working  fluid,  in  this  case  at  the 
heater,  to  increase  the  power  output  due  to  increased  working  fluid  mass  flow  and  also 
by  increasing  the  working  fluid's  specific  heat  The  theoretical  performance  of  different 
variations  of  the  STIG  has  been  evaluated  by  several  authors  [4-17],  It  has  been  known 
for  some  time  that  water/steam  injection  increases  efficiency  and  power  output  [4], 
Renewed  interest  started  during  the  recent  energy  crisis  period  of  a  decade  ago.  Some  of 
the  earliest  work  on  water/steam  injection  was  performed  by  the  NACA  with  the 
objective  of  augmenting  the  performance  of  gas  turbine  aircraft  engines  [5].  Early 
studies  included  water  injection  for  power  increase  only,  which  typically  decreased  cycle 
efficiency  2%  for  each  1%  increase  in  turbine  mass  flow  [6].  The  literature  referenced 
here  studies  various  other  configurations  of  steam  injection  (5,  7-11]  to  include  steam 
injection  after  the  combustor  [7],  water  injection  by  an  evaporator  prior  to  preheating  in  a 
heat  exchanger  [8],  and  steam  generated  by  turbine  exhaust  which  is  used  in  both  a  gas 
turbine  cycle  and  a  steam  turbine  cycle  [9]. 

The  mentioned  advantages  of  the  STIG  over  the  GT  include  the  increases  in  cycle 
efficiency  and  power  output.  Boyce  et  al.  [  10]  describe  a  2-3%  increase  in  efficiency 


and  a  12%  increase  in  power  output  for  each  5%  increase  by  weight  injection  of  steam  in 
their  external  combustion  steam  injected  gas  turbine.  Stochl  [4]  finds  increases  in 
efficiency  and  power  output  as  high  as  30%  and  50%,  respectively,  at  a  turbine  inlet 
temperature  of  1366  K  (2000^).  Even  more  dramatic  is  Fraize  and  Kinney's  [5] 
findings  of  efficiency  increases  of  8  to  14  points  while  increases  for  power  output  vary 
as  much  as  70%  to  130%  for  their  coal-fired  gas  turbine  power  cycle  with  steam 
injection.  This  also  demonstrates  that  STIG  applications  may  be  possible  for  a  variety  of 
fuels  and  not  just  "clean"  fuels.  STIG  systems  have  the  advantage  of  simplicity  when 
compared  to  the  IGT  or  conventional  combined  cycles  of  gas  and  steam  turbines. 
Existing  gas  turbine  cycles  have  been  modified  for  steam  injection  to  handle  higher  loads 
with  no  increase  in  rotating  machinery  [12].  The  STIG’s  simplicity  also  means  smaller 
capital  costs  in  comparison  to  more  complex  systems  Steam  injection  is  known  to 
reduce  the  level  of  pollutants  when  fuel  is  combusted.  Field  tests  have  found  that  a  5% 
by  weight  injection  of  steam  will  reduce  the  amount  of  NOx  emissions  to  acceptable 
levels  [10].  The  STIG  cycle  has  the  disadvantages  of  corrosion  and  water  consumption 
which  all  steam-based  systems  must  cope  with.  Currently,  two  U.S.  commercial  firms 
offer  STIG  cogeneration  systems.  They  are  International  Power  Technology  (TPT),  Palo 
Alto,  California  and  Mechanical  Technology,  Inc.  (MTI),  Latham,  New  York  [13].  IPT 
has  two  operational  installations  of  their  Cheng-Cycle  cogeneration  STIG  systems  at  San 
Jose  State  University  and  a  Sunkist  Growers  installation  in  Ontario,  California  [13,14], 
There  are  no  utilities  currently  operating  steam  injected  gas  turbines,  but  Pacific  Gas  and 
Electric  Co.  is  interested  in  seeing  the  technology  developed  [  13].  Actual  testing  of  a 
STIG  system  in  an  electrical  power  generation  configuration,  with  steam  injected  into  a 
Westinghouse  191-6  gas  turbine  (compression  ratio  at  6.5,  steam/air  mass  ratio  at  .05). 
yielded  a  20%  increase  of  generated  electricity  [15].  After  3000  hours  of  operation, 
there  was  no  ill  effect  on  the  turbine  blades. 


This  study  will  vary  turbine  inlet  temperature  and  steam/air  mass  ratio  (for  the 
STIG)  for  constant  values  of  compression  pressure  ratio  to  evaluate  their  effect  on  cycle 
efficiency  and  power  output.  The  compression  ratios  (denoted  as  Pr)  investigated  are  4, 
8,  12,  16,  20.  This  range  is  common  to  the  studies  in  the  literature.  Brown  and  Cohn 
[16]  state  that  the  highest  pressure  ratio  currently  foreseen  for  single-shaft  industrial  gas 
turbines  is  16.  Although  a  pressure  ratio  of  20  would  provide  the  highest  efficiency,  a 
pressure  ratio  of  16  was  selected  as  the  superior  economic  candidate  for  the  basic  STIG 
cycle.  The  turbine  inlet  temperature  range  investigated  is  1000-2500  K  (1340-4040°F). 
The  temperatures  examined  in  the  literature  were  commonly  in  the  range  of  1073-1700 
K  (1471-2600°F)  with  Day  and  Kidd  [6]  examining  the  range  of  1228-2450  K 
(1750-3950°F).  The  investigated  range  is  admitted  to  be  high.  The  high  end  of  which 
may  not  be  technically  possible  without  exotic  materials  and  elaborate  turbine  cooling 
schemes,  but  an  examination  of  the  effect  of  high  temperatures  on  performance  is  the 
purpose  of  this  study.  The  range  of  steam/mass  air  ratios  (denoted  as  Mr)  investigated 
is  0  to  0.5.  It  will  be  shown  here  that  maximum  First  Law  efficiency  is  achieved  within 
this  range.  This  is  confirmed  by  the  literature. 

The  Second  Law  efficiency  will  also  be  investigated  for  the  configurations  shown 
in  Figure  1  plus  two  GT  cogeneration  applications  where  all  the  steam  produced  by 
turbine  exhaust  goes  to  process  applications  and  is  not  injected  in  the  gas  turbine  cycle. 
Second  Law  efficiency  will  be  defined  here  as  the  ratio  of  the  net  work  output 
determined  by  the  First  Law  analysis  and  the  maximum  reversible  work  determined  with 
the  thermodynamic  states  calculated  in  the  First  Law  analysis.  This  definition  draws  on 
the  work  of  Boehm  [18],  Moran  [19],  and  Kotas  [20],  Each  of  the  systems  to  be 
investigated  here,  the  GT,  IGT,  and  STIG  cycles,  do  not  have  combustion  processes. 
These  systems  have  heat  addition  processes  from  an  external  heat  source.  This  source 
can  be  the  solar  central  receiver  and  its  initial  working  fluid.  This  study  utilizes 


definitions  from  Boehm  [18]  and  Moran  [19]  to  express  the  heat  addition  processes  over 
varying  temperatures  as  reversible  Carnot  equivalent  work.  In  addition  it  is  necessary  to 
complete  the  expression  for  reversible  work  by  adding  the  availabilities  of  the  inlet 
streams  and  subtracting  the  availabilities  of  the  exit  streams.  The  work  of  all  three  of 
these  authors  was  referred  to  in  expressing  these  availabilities. 


CHAPTER  n 


THEORY 

Computer  modelling  of  a  First  Law  analysis  is  the  basis  of  this  study.  The  model 
determines  the  thermodynamic  states  at  each  point  in  the  three  configurations.  (The 
FORTRAN  program  GT  developed  for  this  study  is  in  Appendix  A).  Given  all  the 
necessary  states,  the  program  calculates  the  performance  characteristics  of  each 
configuration:  specific  work  output,  First  Law  efficiency,  and  Second  Law  efficiency. 
This  chapter  will  detail  the  modelling  techiques  used  in  the  computer  program  to 
determine  the  end  states  of  various  processes  and  phenomena.  The  assumed  values  used 
in  the  program  are  listed  in  Table  1. 

Compressor/T  urbine 

The  compression  or  expansion  ratios  and  the  assumed  values  of  the  component 
isentropic  efficiencies  are  the  basis  of  the  calculation  of  the  process  (temperature  change) 
through  these  components.  They  are  assumed  adiabatic.  The  assumed  isentropic 
efficiencies  are  on  the  high  side  of  those  listed  in  the  literature.  For  the  IGT,  the  two 
stages  of  compression  and  expansion  have  the  same  pressure  ratio  which  was 
determined  by  the  square  root  of  the  total  pressure  ratio.  The  final  turbine  expands  the 
working  fluid  to  atmospheric  pressure.  The  calculation  for  end  temperature  is  shown 
below. 

Compm  ssor 

T2s  =  ti  (Pr)k'1/k 


(1) 


Parameter 


Assumed  Value 


Ambient  temperature,  T  j 

300  K 

Ambient  pressure, 

101.325  kPa 

Air  specific  heat  (initial  iteration  value),  cpaQ 

1.0035  kJ/kg  K  [21] 

Air  ratio  of  specific  heats  (initial  iteration  value),  ^ 

1.397  [9] 

Air  molecular  weight,  MWa 

28.97  kg/kmol  [21] 

Air  ideal  gas  constant,  Ra 

.28700  kJ/kg  K  [21] 

Steam  specific  heat  (initial  iteration  value),  cps  ^ 

1.8723  kJ/kg  k  [21] 

Steam  ratio  of  specific  heats  (initial  iteration  value), 

1.327  [21] 

Steam  molecular  weight,  MW§ 

18.015  kg/kmol  [21] 

Steam  ideal  gas  constant,  Rs 

.46152  kJ/kg  K  [21] 

Water  density,  PHjO 

977.0  kg/m3  [24] 

Universal  gas  constant,  R 

8.31434  kJ/kmol  K  [21] 

Compressor  efficiency,  t]c 

0.87 

Turbine  efficiency,  Tlt 

0.89 

Pump  efficiency,  rip 

0.70 

Compressor  pressure  ratios,  Pr 

4,  8,  12,  16,  20 

Pressure  lost  constant,  kpjost 

12.128  (kPa)2/K 

T 


Ti+T2,-T. 

1 2actual  ~  - 

Me 

(2) 

Turbine 

T2s  =  Ti/(Pr)k‘1/k 

(3) 

factual  =  ^1  '  ^t  (^1  ■  ^2s)- 

(4) 

The  specific  work  required  and  produced  by  the  compressor  and  turbine 
processes,  respectively,  is  calculated  as 

w  =  cp  factual '  ^l)  •  (5) 

The  program  does  the  calculations  in  SI  units.  The  specific  work  is  in  terms  of 
kilojoules  per  kilogram  of  compressor  air  flow.  The  calculation  of  the  turbine  specific 
work  as  shown  above  is  additionally  multiplied  by  (0.95  +  Mr).  The  0.95  designates  the 
95%  of  the  air  flow  not  lost  from  the  working  fluid  for  the  turbine  cooling  (more  on  this 
topic  later).  The  steam/air  mass  ratio,  Mr,  is  added  to  the  working  fluid  for  the  ST1G 
configuration.  A  variable  specific  heat  for  air  or  a  steam/air  mixure  working  fluid  is 
determined  using  a  variable  air  specific  heat  relationship  [22]  and  a  variable  steam 
specific  heat  relationship  [23].  These  relationships  use  an  average  temperature  of  the 
start  and  end  points.  For  the  STIG  cycle  the  mixture  specific  heat  is  determined  using 
the  steam/air  mass  ratio,  Mr. 


Heat  and  Reheat 

The  turbine  inlet  temperature  is  an  independent  parameter  in  this  study  so  the  end 
temperature  of  the  heating  process  is  known.  The  specific  heat  input  required  is 
calculated  in  similar  fashion  as  the  work  calculation  using  an  energy  balance  as  shown 


below. 


%  ~  +  Mr)  Cp  (Ttur{jjne  jniet  -  Tj) 


(6) 


Heat  Exchangers 

The  heat  exchangers  to  be  modelled  in  this  study  are  the  intercooler,  the 
regenerator,  and  the  heat  recovery  steam  generator  (HRSG).  A  simple  but  appropriate 
way  was  sought  to  express  their  efficiencies  without  setting  any  values  for  the  heat 
exchanger's  particular  dimensions.  A  similar,  analogous  relationship  was  determined 
and  applied  to  each  heat  exchanger.  This  was  accomplished  by  expressing  heat 
exchanger  effectiveness  in  terms  of  the  ratio  of  the  minimum  temperature  difference 
between  the  two  streams  and  the  total  convective  heat  resistance  of  the  two  streams.  The 
conductive  heat  resistance  of  the  heat  exchanger  is  neglected.  This  ratio  for  each  heat 
exchanger  is  set  equal  to  the  ratio  of  one  representative  heat  exchanger.  This  is  shown 
below. 


AT, 


mm 


R, 


AT, 


nun 


each 
I  HEX 


representative 

HEX 


(7) 


This  can  be  further  expressed  as 


AT 


min 


=  R, 
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HEX 


min 
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HEX 


representative 

HEX 
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— L  +  _JL 

hiAi  h2Aj 


=  AT, 


0110  1  +  1 


h3A2  h4A2 


For  an  assumption  of  equal  heat  transfer  areas,  the  equation  becomes 


representative 


-L  +  _L 

hi  h- 


-L  +  _1 

h-i  hA 


The  representative  minimum  temperature  difference  selected  was  50°C  for  the 
regenerator.  A  value  of  51.1  *C  was  determined  from  a  regenerator  efficiency  of  0.81 
[  16]  for  a  turbine  inlet  temperature  of  1700  K  and  compression  pressure  ratio  of  16.  The 
representative  values  of  the  applicable  convective  heat  transfer  coefficients  (Table  2) 
were  selected  [24,  25],  The  minimum  temperature  difference  between  the  streams  of  the 
other  two  heat  exchangers  was  calculated  and  is  shown  in  Table  2. 

The  locations  of  these  ATmjn  points  are  shown  in  Figure  2.  The  location  is 
determined  by  the  relative  values  of  mass  flow  and  specific  heat  between  the  two 
streams.  These  parameters  determine  the  slope  of  the  two  streams  (shown  as  linear  on  a 
T-length  (x)  diagram)  and  will  determine  where  the  two  streams  are  closest.  The 
locations  of  the  intercooler  and  regenerator  are  pointed  out  but  the  location  of  the  "pinch 
point"  in  the  HRSG  can  be  at  various  locations  due  to  varying  turbine  exhaust 
temperature  and  steam/air  mass  ratio.  This  will  be  discussed  later  in  this  chapter. 


A  method  to  express  pressure  losses  without  using  any  apparatus  parameters  and 
which  was  applicable  at  any  location  was  sought.  For  steady  flow  through  the  com¬ 
ponents,  the  fluid  velocity  and  geometry  of  the  apparatus  is  generally  constant.  There- 
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Medium 

Convective  Heat  Transfer  Coefficient  (W/m^K) 

Forced  convection  of  air 

200 

Forced  convection  of  water 

2000 

Convection  of  boiling  water 

5000 

Heat  Exchanger 

Minimum  Temperature  Difference  (*C) 

Regenerator  (Air- Air) 

50 

Intercooler  (Air-Water) 

27.5 

Regenerator  (Air- Air) 
Intercooler  (Air-Water) 
HRSG  (Air-Boiling  Water) 


26.0 


fore,  the  Reynold's  number  is  constant. 


Re  = 


pVd 

- —  =  constant  . 


(10) 


For  a  constant  Reynold's  number  and  dimensions,  the  friction  factor  will  be  constant. 
The  relationship  for  friction  of  internal  turbulent  flow  is  modified  using  an  expression 
for  mass  flow  rate  and  constant  area.  The  following  equation  for  mass  flow  rate  is  used. 


m  =  pAV  .  (11) 

The  equation  for  the  friction  factor  of  internal  turbulent  flow  is  expressed  as 

friction  =  f(Re)  =  =  constant.  (12) 

pV2 

Applying  Equation  (11),  Equation  (12)  becomes 


App 


=  Ci  . 


nr 


The  ideal  gas  relation  is  used. 


_P 

Ap  RT 


=  C,  . 


nT 


Since  m  and  R  are  constant,  the  equation  can  be  expressed  as 


(Ap)p 

T 


=  C2  . 


(13) 


(14) 


(15) 


This  yields  an  expression  for  pressure  drop  in  terms  of  pressure,  temperature,  and  a 


constant.  A  representative  value  of  a  pressure  drop  of  5  psi  at  atmospheric  conditions 
was  used  to  determine  the  constant,  C2-  The  pressure  drop  through  the  intercooler, 
regenerator,  heaters,  and  HRSG  can  now  be  determined  using  entry  pressure  and  an 
average  temperature. 

The  literature  used  a  number  of  techniques  to  model  pressure  losses  in  one  or 
more  of  the  major  heat  transfer  components.  Fraize  and  Kinney  [5]  set  the  pressure  drop 
across  both  the  combustor  and  steam  generator  at  2%  of  the  initial  pressure.  Bhutani  et 
al.  [9]  applied  a  pressure  drop  of  4%  across  only  the  heater.  Boyle  [11]  accounted  for  a 
4%  pressure  drop  across  the  combustor,  a  4%  drop  in  the  turbine  exhaust  steam  in  the 
HRSG,  and  a  12%  drop  in  the  steam  flow  in  the  HRSG.  He  assigned  generally  larger 
values  to  account  for  parasitic  pipe  losses  also. 

The  method  selected  for  this  study  to  account  for  pressure  losses  has  the 
advantages  of  determining  the  losses  by  use  of  flow  parameters  alone.  This  method  is 
applicable  to  any  of  the  mentioned  components. 

Heat  Recovery  Steam  Generator  (HRSG) 

The  effect  of  varying  steam/air  mass  ratio  and  turbine  exhaust  temperature  on  the 
performance  of  the  HRSG  is  constrained  by  a  minimum  pinch  point  limit  (ATpp) 
between  the  two  streams.  The  amount  of  heat  to  be  recovered  effectively  by  the 
incoming  water  requires  that  the  temperatures  of  the  exhaust  gas  stream  must  be  higher 
than  the  water/steam  stream  temperature  by  a  minimum  amount  at  all  points  in  the  heat 
exchanger.  Another  temperature  difference  to  be  defined  is  the  approach  temperature 
difference  (ATap).  This  is  the  temperature  difference  between  the  entering  hot  exhaust 
stream  and  the  exiting  cooler  water/steam  stream.  The  location  of  the  approach 
temperature  difference  as  well  as  the  possible  locations  of  the  pinch  point  temperature 
difference  are  shown  in  Figure  3.  The  literature  describes  the  efficiency  of  the  HRSG  by 
setting  ATap  at  a  certain  value  and  then  checks  to  see  if  the  ATpp  minimum  limit  has  not 


Turbine  inlet  temperature 
is  not  hot  enought  to  boil 
water  in  HRSG. 


Turbine  inlet  temperature 
is  hot  enough  to  start  water 
to  boil.  Steam  exits  HRSG 
in  saturation. 


Turbine  inlet  temperature 
is  high  enough  for  steam 
to  reach  superheat. 


re  3.  Possible  heat  recovery  steam  generator  (HRSG)  pinch  point  (pp)  locations. 
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been  violated  at  the  point  where  boiling  of  the  water  begins  (point  6  in  Figure  3).  This 
sets  the  value  of  the  exiting  steam  to  the  temperature  of  the  turbine  exhaust  minus  ATap 
Once  the  temperature  difference  at  the  beginning  of  boiling  approaches  ATpp,  the 
temperature  of  the  exhaust  stream  at  the  pinch  point  is  set  at  the  water's  saturation 
temperature  plus  ATpp.  The  temperature  of  exiting  steam  is  then  allowed  to  vary. 

The  range  of  values  of  ATpp  in  the  literature  is  16.7  -  27.8*C  (30  -  50*F)  [4,  5,  7, 
9, 11,21,  26].  Stochl  [4]  considers  the  value  of  27.8’C  (50*F)  as  a  reasonable  value  for 
both  ATap  and  ATpp.  In  his  study  he  examined  four  values  for  each.  Other  authors 
have  set  ATap  at  a  different  value  than  ATpp  including  41.7‘C  (75*F)  [16]  and  55.5'C 
(100*F)  [9].  The  value  of  26*  C  developed  in  this  study  for  ATpp  is  certainly  in  line  with 
the  literature.  In  terms  of  the  literature  this  study  places  both  ATap  and  ATpp  at  26‘C. 
Figure  3  shows  that  the  steam  exiting  the  HRSG  can  be  superheated,  saturated,  or  can 
still  be  liquid  depending  on  the  turbine  exhaust  temperature  and  the  steam/air  mass  ratio. 
The  modelling  of  the  HRSG  in  this  study  can  handle  each  possibility.  The  modelling 
done  in  the  literature  commonly  describes  the  performance  of  the  HRSG  only  when 
^ap  is  set  for  exiting  superheated  steam  and  when  the  ATpp  value  is  applied  at  the  first 
point  of  boiling  for  steam  exiting  as  superheated  steam  and  saturated  steam.  In  Figure  3 
these  three  situations  are  designated  as  pinch  points  5,  6,  and  3,  respectively.  It  will  be 
shown  in  the  results  of  this  study  that  the  steam/air  mass  ratio  which  achieves  the 
maximum  First  Law  efficiency  is  the  ratio  which  causes  the  pinch  point  (minimum 
temperature  difference  between  the  two  streams)  to  move  from  the  approach  temperature 
location  to  the  first  point  of  boiling.  This  point  of  maximum  First  Law  efficiency  is 
confirmed  by  the  literature.  To  describe  the  STIG  performance  over  the  entire  range  of 
turbine  inlet  temperatures  and  steam/air  mass  ratios  the  modelling  of  the  HRSG  was 
developed  to  handle  all  seven  possible  pinch  point  locations. 

The  known  conditions  within  the  HRSG  before  calculations  are  made  include  the 
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turbine  exhaust  temperature,  the  ratio  of  steam  to  air  in  the  turbine  exhaust,  inlet  water 
temperature,  and  the  water  saturation  temperature.  The  water  is  pumped  to  compressor 
discharge  pressure  plus  the  calculated  pressure  loss  through  the  HRSG  and  combustor. 
The  saturation  temperature  for  that  pressure  at  each  of  the  selected  compression  ratios 
and  the  other  saturation  properties  of  latent  heat  of  vaporation  (hfg)  and  Sfg  (for  the 
Second  Law  analysis)  were  interpolated  from  Keenan  et  al.  [27].  Additionally,  a  liquid 
water  specific  heat  value  for  each  compression  ratio  was  determined  by  selecting  the 
corresponding  value  of  Cp  for  an  average  temperature  (between  ambient  and  saturation 
conditions)  from  the  curve  in  Figure  2  in  Keenan  et  al.  [27].  The  two  unknown 
temperatures  are  the  exiting  steam  temperature  and  the  exiting  hot  stream  exhaust 
temperature  (denoted  as  Tstac^).  These  unknown  temperatures  are  found  by  an  iterative 
process.  First  the  pinch  point  temperature  difference  (ATpp)  of  26°C  is  set  at  first  two 
possible  locations  (points  1  and  2  in  Figure  3)  and  the  unknown  temperatures  are  found 
by  equating  the  heat  given  up  by  the  hot  stream  to  the  heat  added  to  the  cold  stream. 
This  assumes  no  heat  loss  from  the  HRSG.  In  both  cases  the  temperature  difference  at 
the  other  possible  pinch  point  location  is  checked  to  see  if  it  is  more  applicable.  If  the 
heat  transferred  is  sufficient  to  heat  the  water,  for  that  stream/air  mass  ratio,  to  saturation 
then  the  next  two  possible  pinch  point  locations  (points  3  and  4  in  Figure  3)  for  exiting 
saturated  steam  are  checked.  If  the  heat  transferred  in  both  of  those  situations  is  enough 
to  heat  the  steam  to  superheat  then  the  last  three  possible  pinch  point  locations  for  this 
case  (points  5, 6,  and  7,  Figure  3)  are  investigated  to  see  which  is  applicable.  In  general 
the  specific  heat  (kJ/kg  air)  given  up  by  the  hot  stream,  which  is  equated  to  the  total  heat 
transferred,  is  calculated  as  shown  below. 

^HEX  =  (0.95  +  Mr)  •  Cp^^^e  •  (Turbine  exit  *  ^stack  )  (^) 

The  specific  heat  added  to  the  cold  stream  is  determined  as  shown  below  (in  the  case  of 
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superheat  steam  exiting). 


q  -  Mr  [cpWater  •  (Tsat  -  Tj)  +  hfg  +  cRteam  •  ^exiting  -  Tsat)]  . 

steam 


(17) 


The  HRSG  mcxielling  subroutine  also  finds  the  required  heat  input  to  raise  the  exiting 
steam  to  the  desired  turbine  inlet  temperature.  This  quantity  is  found  by  subtracting  the 
heat  transferred  in  the  HRSG  from  the  total  heat  required  to  raise  the  water  from  ambient 
conditions  to  the  turbine  inlet  temperature.  This  is  shown  below. 

q|Hteam*  Mr  ,CPwater  <Tsa< '  Tl)  +  hfg  +  cPsteam'  turbine '  Tsat»  ’  <)HEX 

Met  (18, 


As  mentioned  earlier  the  values  of  specific  heat  at  constant  pressure  (Cp)  for  the  air, 
steam,  and  steam/air  mixture  streams  are  determined  at  average  temperatures  from  tem¬ 
perature  dependent  relations  for  air  and  steam  specific  heat  (Cp  (T))  [22,  23]. 

Pump 

The  required  work  to  pump  the  water  to  the  desired  compressor  discharge  pres¬ 
sure  is  calculated  as 


Specific  Work  = 


Mr  •  (Pr-Pj  +  APhrSG) 
Pwater  •  ^p 


(19) 


This  quantity  is  negligible  for  the  range  of  parameters  of  this  study. 

Turbine  Cooling 

For  higher  temperature  applications  in  gas  turbines  a  certain  amount  of  compressor 
air  is  bled  off  from  the  working  stream  to  cool  the  turbine  blades  to  keep  them  within 
their  materials'  operating  temperature  limits.  There  are  several  modelling  techniques 
listed  in  the  literature.  Fraize  and  Kinney  [5]  state  that  the  amount  of  air  bled  from  the 


compressor  is  1.5%  of  the  total  compressor  flow  for  each  10CFF  above  a  turbine  inlet 
temperature  of  1750°F  (1228  K).  Day  and  Kidd  16]  set  a  schedule  of  4%  flow  loss  for  a 
turbine  inlet  temperature  of  1500°F  (1089  K),  a  8.5%  loss  for  a  temperature  of  I750T 
(1228  K),  and  a  16%  loss  for  a  temperature  of  200CTF  (1367  K).  For  this  study  the 
technique  of  Bhutani  et  al.  [9]  was  used.  This  is  a  straight  5%  air  flow  loss  to  turbine 
cooling.  This  was  reflected  in  earlier  listed  equations.  This  is  an  admittedly  simple 
calculation.  The  authors  of  this  technique  applied  it  to  their  study  which  had  a  1700  K 
(2600°F)  turbine  inlet  temperatue.  They  make  a  good  point  that  in  reality  the  initial 
quantity  of  bleed  flow  would  be  larger  but  certain  portions  of  the  air  coolant  are  regained 
to  the  working  fluid  at  all  subsequent  turbine  blade  stages  to  the  stage  at  which  that  air 
stream  is  applied  for  cooling.  In  keeping  with  the  desire  not  to  be  locked  in  specifying 
any  particular  component  parameters  or  dimensions  than  is  necessary  this  technique  of 
determining  turbine  cooling  bleed  off  was  used. 


The  calculation  of  the  First  Law  efficiency  is  a  straightforward  calculation  of  the 
ratio  of  net  work  output  to  heat  input.  For  the  IGT  cycle  there  are  two  stages  of 
compression  and  expansion  as  well  as  two  heating  components.  The  STIG  cycle  has  a 
work  requirement  for  the  water  pump  and  a  heat  input  requirement  to  heat  the  steam 
exiting  the  HRSG  to  the  turbine  inlet  temperature.  The  First  Law  efficiency  calculation 
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Second  Law  Efficiency 

At  present  there  are  no  universally  accepted  definitions  for  Second  Law  efficiency. 
This  study  will  follow  the  concepts  of  Boehm  [18],  Moran  [19],  and  Kotas  [20]. 
Second  Law  efficiency  will  be  defined  here  as  the  ratio  of  the  net  work  and  the  reversible 
work.  The  reversible  work  is  further  defined  as  the  sum  of  the  inlet  stream  availabilities 
and  Carnot  equivalent  work  of  each  of  heat  inputs  minus  the  sum  of  the  exit  stream 
availabilities. 


w 


rev  -  ;^(b  +  w„-Ih 


inlets 


heat  e  exits 
inputs 


(21) 


The  results  to  be  shown  later  will  show  the  reversible  work  with  and  without  accounting 
for  the  exit  stream  availabilities  as  a  negative  term  in  the  reversible  work.  It  may  be 
argued  that  the  exit  streams  that  are  vented  to  the  stack  are  not  available  and  need  not  be 
assessed  as  a  loss  or  penalty.  Some  further  definitions  need  to  be  made.  First, 
thermodynamic  availability  is  defined  as  the  amount  of  work  that  can  be  produced  when 
a  substance  at  a  given  state  reverts  to  the  "dead"  state  while  producing  work  with  a 
Carnot  cycle  [18].  Here  the  dead  state  is  usually  defined  as  atmorpheric  conditions  and 
the  hypothetical  Carnot  cycle  is  assumed  to  be  rejecting  heat  to  the  ambient  temperature. 
T  j.  For  this  study  the  availability  of  a  stream  is  expressed  as  showm  below. 


b  =  h-h1  -Tj  (s  -  sj)  .  (22) 

The  kinetic,  chemical,  and  potential  components  of  the  availability  are  considered  negli¬ 
gible.  For  a  mixture  of  n  ideal  gases,  the  availability  becomes  (in  molal  basis)  [19,  20] 


rT  Ti 

b=k=l  yk  [  1  (1‘T)CPk 

T1 


dT]  +  R  T  j  In  (— ) 
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The  Carnot  equivalent  work  of  the  heat  inputs  is  calculated  as 


T1 

w  =  (l-T)qh  • 


(24) 


In  comparison  to  this  study's  definition  of  Second  Law  efficiency,  the  definitions  pro¬ 
vided  in  the  literature  are  provided  here.  Moran  [19]  provides  an  example  of  a  heat 
engine  where  work  is  produced  by  drawing  heat  from  a  high  temperature  reservoir,  at 
temperature  Tpj,  and  rejecting  heat  to  a  low  temperature  reservoir,  at  temperature  T^. 
The  Second  Law  efficiency  can  be  expressed  as  follows  (TQ  is  ambient  temperature): 
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e  = 


net 
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He  states  that  this  expression  does  not  charge  as  a  loss  the  availability  exiting  in  the  heat 
interaction  at  T^.  This  study  follows  more  closely  this  alternate  definition,  which  as 
he  describes  charges  the  availability  exiting  in  the  heat  interaction  at  T^  as  a  loss.  This  is 
expressed  as 


W 


e  = 


net 


(I-^SQh 


(26) 
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This  study  used  this  approach  and  assesses  availabilities  entering  and  leaving  the  con¬ 
trol  volume  of  the  system.  Kotas  [20]  provides  a  Second  Law  efficiency  definition  for  a 
gas  turbine  plant.  This  plant  bums  fuel  internal  to  the  cycle.  His  Second  Law  efficiency 
is  expressed  below. 
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e  _ _ Rated  plant  net  output _  2 

Total  availability  of  all  input  fuel  components 

Kotas  and  El-Masri  [28]  provide  details  of  the  use  of  finding  stream  availability  at  each 
of  the  cycle  state  points  to  locate  the  cause  and  magnitude  of  losses  within  the  system. 
This  study  uses  a  control  volume  approach  to  evaluate  overall  cycle  performance,  so 
finding  the  location  and  magnitude  of  the  system's  irreversibilities  is  not  a  goal  of  this 
study.  This  approach  of  a  Second  Law  analysis  is  still  quite  new,  but  has  a  promising 
future. 

Simple  Bravton  Cycle  (GT) 

For  the  GT  and  IGT  cycles  the  referenced  dead  state  is  air  at  ambient  temperature 
and  pressure.  The  availability  of  the  air  inlet  stream  is  then  zero.  The  Second  Law 
efficiency  (e)  for  the  GT  cycle  is  defined  as  shown. 


wnet 

weair  ~bexhaust 


(28) 


The  Carnot  equivalent  work  of  the  heat  input  is  calculated  as  shown  below. 


(•t3*  T!  [T3  T, 

wCajr=0.95j  (1  -— )  Sq -0.95  J  (1  -  -)cpairdT  (29) 

T 

W9aii  =  0  95  >ir  I(T3  '  T2> '  T1  ln  tjr  »  ■  <30> 

The  exhaust  availability  is  calculated  as 

bexhaust  =  °-95[Ah-Tl(As)]-  <3]) 

♦Refer  to  Figure  1  to  find  the  locations  of  referenced  state  points. 
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Since  P4  =  Pj,  then  this  term  can  be  expressed  as 


3exhaust  ~  ^-95  cp^j.  '  ^i)  -  Tj  in  ( ^  )] 


As  mentioned  above,  the  availability  of  the  inlet  air  stream  is  zero.  In  the  IGT 
cycle  there  are  two  heating  processes  and  the  cooling  of  the  air  stream  in  the  intercooler. 
These  are  expressed  as  three  Carnot  equivalent  heat  inputs.  Each  is  calculated  as  shown 
above  for  the  GT  cycle.  The  exhaust  gas  availability  is  determined  using  the  exhaust 
temperature  from  the  regenerator.  The  cycle  Second  Law  efficiency  is  defined  as 


i?l  Weheatingj  We  intercooling  Exhaust 


A  Second  Law  efficiency  for  the  STIG  cycle  is  computed  in  an  analogous  manner 
to  the  approach  covered  for  the  GT  and  IGT  cycles.  In  addition,  two  other  Second  Law 
efficiencies  will  be  developed  for  comparison  in  which  all  the  steam  produced  in  the 
HRSG  goes  to  process  work  elsewhere  and  is  not  injected  into  the  gas  turbine.  In  these 
second  two  cases,  one  will  heat  the  exiting  steam  to  the  turbine  inlet  temperature  and  the 
other  will  not  heat  the  steam  beyond  its  HRSG  exit  temperature. 

In  evaluating  the  Second  Law  efficiencies  for  these  three  systems,  it  will  allow  a 
comparison  of  a  steam  injectt  i  system  and  a  cogeneration  system.  These  systems  are 
alike  except  the  STIG  system  injects  the  steam  produced.  An  evaluation  of  the  level  of 


irreversibilities  of  the  steam  injection  can  be  made. 

For  the  first  definition  of  the  Second  Law  efficiency,  the  reference  state  is  the 
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steam/air  mixture  at  ambient  temperature  and  pressure.  In  this  case,  the  inlet  air  and 
water  streams  are  at  ambient  pressure  which  is  greater  than  their  corresponding  partial 
pressure  in  the  exhaust  mixture  stream.  The  inlet  stream  availabilities  are  then  nonzero. 
Using  the  earlier  definition  of  availability  (Equation  (23)),  the  availability  of  the  air  inlet 
stream  (on  a  mass  basis)  can  be  expressed  as 

Tj  *  R  In  (1/y  air) 

a^r  inlet  MWair 

The  water  inlet  stream  is  expressed  in  a  similar  fashion.  The  availability  of  the  mixture 
exhaust  is  expressed  as 

Tctsrk 

^xhaus.  =  (Mr  +  °'95)  •  Cftnix  ‘  [T«**  -  T,  -  T,  •  In  )].  (35) 

The  heat  input  requirements  for  the  STIG  cycle  are  to  the  air  and  the  steam.  Depending 
on  the  turbine  exhaust  temperature  and  the  steam/air  mass  ratio,  the  steam  can  exit  the 
HRSG  at  superheat  or  in  saturation.  The  Carnot  work  equivalent  for  the  air  heat  input  is 


The  equivalent  work  for  the  heat  input  to  take  the  steam  through  superheat  to  the  desired 
turbine  inlet  temperature  is  calculated  as 


"‘superheat  = tb' '  CPsteam  [  ^  '  T*>') '' T>  'n  (^t  »  ' 
steam 


(39) 


The  Second  Law  efficiency  for  this  application  is  defined  as 
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(40) 


In  the  last  two  Second  Law  efficiencies  to  be  described  for  the  GT  cycle  with  stream 
production,  the  air  and  water  streams  never  mix.  The  reference  states  are  separate  air 
and  water  streams  at  ambient  conditions.  The  availability  of  these  inlet  streams  is  zero. 
The  exit  stream  of  steam  is  produced  to  do  process  work.  Its  availability  will  be  added 
to  wnet  in  the  numerator  of  the  efficiency  ratio.  For  this  cycle,  it  was  found  that  the 
steam  could  exit  as  liquid  water,  saturated  steam,  or  superheated  steam.  This  is  due  to 
the  loss  of  increased  turbine  exhaust  specific  heat  when  steam  is  not  injected  into  the 
working  fluid.  With  only  air  in  the  hot  stream  of  the  HRSG,  the  HRSG  modelling 
subroutine  must  be  able  to  handle  all  7  possible  pinch  points  as  described  earlier.  The 
availability  of  the  HRSG  exhaust  stream  is  done  as  shown  before  except  it  is  for  only 
air.  If  the  steam  exits  as  liquid  water  its  availability  is  described  as 


1  waterexit 

^water =  cRvater  ^ water  *  Tj)  -  Tj  In  (  ~  )] 

exit 


1 


(41) 


If  the  steam  exists  as  saturated  steam,  the  quality  of  the  steam  is  found  to  properly  ac- 


count  for  the  (h  -  hf)  and  (s  -  Sf)  contributions  to  the  availability  by  expressing  them  as 
x  •  hfg  and  x  -  Sfg.  The  quality  is  found  as 

h-hf  (h  -  hj)  -  (hf  -  h1)  qnEX  *  ^sat 


The  availability  is  then  expressed  as 


^sat  steam  ^  {  cpwater  t^sat '  ^i)  *  "^1  (  y  )1  +  x  1  hfg  *  Tj  •  x  •  Sfg}. 


For  steam  that  exits  as  superheated  steam  (and  similarly  the  case  where  the  steam  is 
heated  to  the  turbine  inlet  temperature)  the  availability  becomes  the  following  relation. 


1  Sclt 

^superheat =  ^cPwater  l^sat  *  Tj)  *  Tj  In  (  ~ y  )]  +  hfg  -  Tj  Sfg 
steam  1 


+  cPsteam  ‘  l^exit  *  ^sat^  '  In  (  M  •  (^4) 

*sat 

Additionally  for  the  case  where  the  steam  is  heated  to  turbine  inlet  temperature,  the  cal¬ 
culation  of  equivalent  work  of  heat  input  must  also  handle  the  three  possibilities  of  exit¬ 
ing  steam  (water,  saturated  steam,  and  superheated  steam).  The  calculation  for  the  case 
when  the  steam  exits  as  liquid  water  will  show  how  each  case  is  done. 
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The  Second  Law  efficiency  for  the  case  of  no  steam  injection  and  where  the  steam  is 
heated  to  turbine  inlet  temperature  is  expressed  as 


wnet  +  ^steam(turbine  inlet  temp) 
Weair  +  Westeam '  ^^exhaust 


(46) 


The  Second  Law  efficiency  for  the  case  of  no  steam  injection  and  where  the  steam  is  not 
heated  beyond  HRSG  exit  temperature  is  expressed  as 

e  _  wnet +  ^steam(HRSG  exit  temp) 

Weair '  ^rexhaust 

This  method  of  expressing  the  Second  Law  efficiency  for  these  GT  cycles  with  steam 
production  closely  follows  the  Second  Law  efficiency  definition  provided  by  El-Masri 
[28].  His  definition  is 
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(48) 


Additionally  he  provides  an  interesting  expression  of  the  Second  Law  efficiency  in  terms 
of  the  Carnot  efficiency  and  irreversibilities  at  each  stage  of  the  cycle. 


e  ^Carnot  •  ^  L  (49) 

components 


He  expresses  L  as  the  ratio  of  the  stream  availability  loss  in  a  process  to  the  total  heat 
added  to  the  cycle. 

This  study  presents  several  variations  of  the  definition  of  Second  Law  efficiency. 
As  a  summary,  these  definitions  are  recapped  in  Table  3.  Generally,  each  definition  is  a 
ratio  of  net  work  to  reversible  work.  Depending  on  how  the  reference  state  is  defined, 
the  reversible  work  may  include  inlet  stream  availabilities.  It  will  be  shown  in  the  next 


Second  Law  Efficiency  Definitions 


Cycle  Second  Law  Efficiency  (e)a 
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aDcfiniiions  do  not  include  exhaust  stream  availability  which  will  be  shown  to  be  not  useful  to  this  studs 
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chapter  that  including  the  exit  stream  availability  in  the  reversible  work  is  not  useful  to 
this  study.  The  major  components  of  the  reversible  work  are  the  Carnot  equivalent  work 
of  the  heat  added  to  the  air  and,  depending  on  the  cycle,  added  to  the  liquid  water, 
saturated  steam,  and  superheated  steam.  For  the  GT  cycle  with  the  HRSG  to  produce 
steam,  the  availability  of  the  steam  at  either  turbine  inlet  temperature  or  HRSG  exit 
temperature  is  added  to  the  net  work. 


The  limits  to  cycle  performance  are  a  result  of  maximum  limits  to  the  parameters  of 
steam/air  mass  ratio  and  maximum  cycle  temperature.  Both  of  these  parameters  are 
influenced  by  several  variables. 


The  appearance  of  condensation  in  the  exhaust  gas  from  the  HRSG  of  the  ST1G 
cycle  places  a  maximum  limit  on  the  steam/air  mass  ratio.  If  the  stack  temperature  at  the 
HRSG  exit  is  below  the  dew  point  temperature  at  ambient  pressure  then  possibly 
corrosive  condensation  will  occur  within  the  HRSG.  Although  not  the  case  in  this 
study,  when  fuel  containing  sulfur  is  burned,  any  condensation  in  the  exhaust  flow  will 
be  acidic.  The  literature  has  selected  a  number  of  safe  minimum  stack  temperature  limits 
to  preclude  the  formation  of  acidic  condensation.  Stochl  [4]  places  a  minimum  Tstack 
limit  in  the  range  of  250-300°F  (394-417  K)  while  Bhutani  et  al.  [9]  placed  the  minimum 
limit  at  275°F  (408  K)  and  Day  and  Kidd  [6]  selected  310°F  (428  K)  as  the  limit. 
Depending  on  the  degree  of  water  treatment  of  the  feedwater  for  the  STIG  cycle,  the 
condensation  of  the  steam  in  the  HRSG  can  have  a  corrosive  potential  for  long  term 
operation  The  effect  on  overall  cycle  performance  to  preclude  any  condensation  within 
the  STIG  cycle  will  be  shown  later. 


Visible  Plume 

Depending  on  the  ambient  conditions,  a  visible  exhaust  plume  from  condensing 
steam  can  form  at  a  stack  temperature  above  the  prevailing  dewpoint  temperature. 
Although  the  plume  is  only  steam,  a  visible  plume  is  a  sign  of  possible  pollution.  To 
preclude  a  visible  plume,  a  lower  limit  than  d(  >ign  conditions  would  be  placed  on  the 
maximum  steam/air  mass  ratio.  In  areas  where  large  fluctuations  of  ambient  conditions 
occur,  it  may  be  necessary  to  install  sensing  equipment  to  monitor  ambient  temperature 
and  relative  humidity  so  that  the  steam/air  mass  ratio  can  be  adjusted  to  eliminate  a  visible 
plume. 

Compressor  Surge  Limit 

In  fixed-geometry  machines  not  specifically  designed  for  steam  injection  the 
relative  increase  in  mass  flow  through  the  turbine  with  steam  injection  can  result  in  an 
increase  in  compressor  back  pressure  and  therefore  an  increase  in  the  compressor 
pressure  ratio.  As  the  steam  injection  increases,  the  compressor  pressure  ratio  increases 
toward  its  surge  limit.  The  maximum  increase  in  performance  obtainable  with  steam 
injection  then  depends  on  the  surge  margin  available  to  the  particular  compressor  [4j. 
This  study  did  not  place  any  limits  on  steam/air  mass  ratio  as  no  particular  compressor 
surge  limits  were  specified. 

Complete  Oxygen  Combustion 

This  limit  pertains  to  cases  when  fuel  is  burned  in  the  combustor  and  is  included 
here  for  information  only.  Complete  oxygen  combustion  serves  as  a  limit  to  both 
maximum  cycle  temperature  and  steam/air  mass  ratio.  For  a  certain  steam/air  mass  ratio 
there  is  no  longer  any  excess  oxygen  to  bum  additional  fuel  to  heat  an  increase  in  steam 
injection  to  the  desired  temperature.  Further  additions  to  steam  are  not  possible  if  the 
turbine  inlet  temperature  is  not  reduced  [6, 1 1]. 
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Materials'  Temperature  Limit 

As  discussed  earlier,  turbine  cooling  is  necessary  for  higher  temperature 
applications  in  gas  turbine  cycles.  Although  over  a  billion  dollars  has  been  spent  since 
World  War  II  on  metallurgical  research  in  developing  new  alloys  to  extend  jet  engine  and 
gas  turbine  temperature  limits,  there  has  only  been  an  advance  of  approximately  lOO'C 
[29].  The  maximum  permissible  metal  temperature  in  gas  turbine  blades  has  been 
extended  from  ~700  to  ~  80CTC  (~1 292  to  ~1472‘F).  Boyle  [11]  sets  a  maximum 
temperature  limit  at  1089  K  (1500°F).  This  temperature  limit  is  based  on  the  fact  that  in 
40  years  of  research  the  same  type  of  alloys  of  Fe  -  Cr  -  Ni  are  still  primarily  used. 
These  alloys  have  a  melting  point  around  1643  K  (2550‘F)  giving  a  sintering 
temperature  (local  welding  at  75%  melting  temperature)  of  -  1230  K  (~  1750°F).  Fraize 
and  Kinney  [5]  state  that  at  the  time  of  their  writing  DOE  programs  were  ongoing  to 
extend  ceramic  technology  to  applications  for  gas  turbines.  This  will  be  a  difficult  task. 
Ceramics  have  excellent  resistance  to  thermal  stresses  but  are  very  brittle.  The  addition 
of  aluminum  to  ceramic  compounds  is  being  tried  to  improve  ductility  [1].  Successful 
ceramic  compounds  could  possibly  extend  maximum  temperatures  to  1561  K  (2350°F) 
[3],  Application  of  coated  carbon-carbon  materials  may  extend  temperatures  to  1700  K 
(2600°F)  [3],  At  this  point  in  time,  turbine  inlet  temperatures  can  be  extended  as  far  as 
the  cooling  schemes  can  keep  the  turbine  blades  within  the  permissible  temperature 
limits.  Fraas  [29]  states  that  with  air  cooling  there  is  a  point  of  diminishing  returns 
where  efficiency  decreases  due  to  increasing  compressor  air  bleed  off  requirements  to 
cool  the  turbine.  Boyce  [1]  reports  that  this  point  of  diminishing  returns  is  1533  K 
(2500*F).  Fraas  [29]  does  place  hope  in  water  cooling  of  gas  turbine  blades  which  may 
allow  turbine  inlet  temperatures  as  high  as  2273  K  (2600"F).  The  water  cooling  will 
require  pumping  which  in  terms  of  cycle  efficiency  will  be  a  much  smaller  cost  than 
bleed  off  from  the  working  fluid.  It  is  therefore  unlikely  that  the  higher  temperature 


limits  of  this  study  (range  is  1000-2500  K)  will  be  applicable  for  some  time  yet. 


Steam  Dissociation 

Another  limit  to  the  maximum  cycle  temperature  is  the  degree  of  dissociation  of 
steam  at  the  increasing  turbine  inlet  temperatures.  Bhutani  et  al.  [9]  placed  a  limit  of 
1  100'F  (1000°F  as  a  practical  limit)  on  the  superheat  of  steam  in  Rankine  systems  due  to 
the  decomposition  of  the  steam.  Fraas  [29]  describes  the  primary  hazard  is  that  steam 
increasingly  dissociates  at  temperatures  above  ~  315“C  (600°F)  which  allows  the 
hydrogen  to  diffuse  through  the  metal  walls  and  frees  the  oxygen  to  attack  the  metal 
walls.  The  literature  referenced  [4-17]  has  not  placed  any  temperature  or  steam/air  mass 
ratio  limits  for  the  STIG  cycle  due  to  this  phenomenon.  In  some  application  of  gas 
turbines  dissociation  can  be  valuable  from  the  standpoint  of  increasing  the  turbine 
volume  flow  relative  to  the  compressor  volume  flow  [29]. 


CHAPTER  III 


RESULTS  AND  DISCUSSION 

This  chapter  will  describe  the  results  of  the  First  Law  and  Second  Law  Analyses 
of  the  performance  of  the  three  Brayton  cycle  systems. 

First  Law  Analysis 

Figures  4  through  8  show  plots  of  First  Law  cycle  efficiency  versus  turbine  inlet 
temperature  for  the  three  systems  at  each  of  the  compressor  ratios.  Figures  9  through  13 
show  the  plots  of  specific  net  work  versus  turbine  inlet  temperature  for  each  pressure 
ratio.  A  fourth  curve  (dash  line)  is  shown  on  Figures  4  through  13  to  depict  the  effect  of 
limiting  the  steam/air  mass  ratio  so  that  no  condensation  occurs  in  the  STIG  cycle.  Each 
of  the  plots  for  the  STIG  cycle  uses  the  steam/air  mass  ratios  which  achieve  the 
maximum  First  Law  efficiency.  This  includes  the  plots  of  specific  net  work  which 
reflect  the  net  work  at  the  steam/air  mass  ratio  which  achieve  the  maximum  efficiency. 
Each  of  the  plots  does  not  cover  the  entire  turbine  inlet  temperature  of  1000-2500  K. 
The  plots  were  cut  off  just  before  the  temperature  where  a  steam/air  mass  ratio  of  0.500 
no  longer  provided  a  maximum  efficiency  for  the  STIG  cycle.  Performance  curves  for 
the  IGT  cycle  were  not  provided  at  the  smallest  pressure  ratio,  Pr  =  4,  since  at  lower 
turbine  inlet  temperatures,  the  2  turbines  cannot  provide  enough  work  to  run  the  2 
compressors. 

On  comparing  the  3  cycles  in  Figures  4  through  13  it  can  be  seen  that  both  the 
Brayton  derivative  cycles  (IGT  and  STIG)  provide  improvements  in  both  efficiency  and 
net  work  output  over  the  simple  Brayton  cycle  (GT).  The  IGT  cycle  is  superior  than  the 
other  two  cycles  in  efficiency  over  the  entire  temperature  range.  The  STIG  cycle 
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provides  a  considerably  greater  improvement  in  net  work  output  than  the  IGT  cycle  over 
much  of  the  temperature  range.  As  the  pressure  ratio  increases  the  STIG  cycle  achieves 
a  maximum  efficiency  at  lower  steam/air  mass  ratios.  At  the  highest  two  pressure  ratios, 
16  and  20,  there  appears  to  be  a  discontinuity  in  the  specific  net  work  plots  for  the  STIG 
cycle.  In  this  area  of  higher  pressure  ratio  and  lower  turbine  inlet  temperature  the 
maximum  efficiency  is  achieved  at  a  different  HRSG  pinch  point  location  than  the  rest  of 
the  studied  regime.  This  pinch  point  location  occurs  at  higher  steam/air  mass  ratios.  The 
higher  steam/air  mass  ratios  yield  higher  net  work  outputs.  This  phenomenon  will  be 
covered  later. 

The  GT  and  IGT  efficiencies  continue  to  increase  with  temperature  but  do  so  with 
a  much  diminished  return.  The  STIG  cycle  as  mentioned  achieves  a  maximum 
efficiency.  The  performance  of  the  STIG  cycle  is  dictated  by  the  performance  of  its  heat 
recovery  steam  generator  (HRSG).  For  the  temperature  range  studied  here,  as  steam  is 
introduced  in  the  HRSG  at  small  mass  ratios  the  turbine  exhaust  stream  is  hot  enough  to 
heat  the  smaller  quanities  of  steam  relative  to  air,  to  a  superheated  vapor.  The  streams 
are  closest  at  pinch  point  #5  as  shown  in  Figure  3.  The  relative  specific  heat  values  (Cp) 
and  mass  flow  ratios  of  the  two  flow  streams  are  such  that  the  slope  of  hot  stream  as 
seen  in  Figure  3  is  quite  shallow.  Pinch  point  5  is  the  appropriate  location  for  the  small 
mass  ratios  and  the  high  temperatures  covered  in  this  study.  As  the  mass  flow  rate 
increases  the  turbine  exit  temperature  increases  due  to  the  increase  in  the  turbine  working 
fluid's  specific  heat.  The  exit  temperature  of  the  steam  leaving  the  HRSG  also  increases 
since  its  temperature  is  fixed  to  the  increasing  turbine  exit  temperature  by  the  minimum 
pinch  point  difference  (ATpp).  With  the  pinch  point  at  this  location,  the  efficiency 
increases  with  the  mass  ratio.  In  other  words,  the  gain  in  the  turbine  work  is  greater 
than  the  penalty  of  heating  more  steam  to  the  desired  turbine  inlet  temperature. 

The  phenomenon  of  increasing  efficiency  with  steam/air  mass  ratio  does  not  go  on 


indefinitely.  The  specific  heat  of  the  turbine  exhaust  increases  with  die  mass  ratio.  The 
slope  of  the  hot  exhaust  in  the  HRSG,  as  shown  in  Figure  3,  decreases  as  the  steam/air 
mass  ratios  increases.  The  difference  in  the  temperature  of  the  two  streams  finally 
reaches  the  ATpp  value  at  the  first  point  of  boiling.  This  location  was  noted  as  pinch 
point  location  #6  in  Figure  3.  With  further  increases  in  the  steam/air  mass  ratio  the 
minimum  temperature  difference  remains  fixed  at  this  location  and  the  exiting  steam 
temperature  will  then  decrease  as  the  mass  ratio  increases.  From  now  on  increases  in  the 
turbine  work  are  less  than  the  penalty  of  heating  increasing  amounts  of  steam  to  the 
desired  turbine  inlet  temperature.  After  the  pinch  point  location  transition  occurs  the 
cycle  efficiency  decreases  as  mass  ratio  increases.  Throughout  much  of  the  studied 
regime  the  mass  ratio,  which  causes  this  transition  to  occur,  provides  the  maximum 
efficiency  plotted  in  Figures  4  through  8.  This  is  verified  by  similiar  plots  in  the 
literature  [4,  6,  9,  19,  11,  16,  26], 

If  the  steam/air  mass  ratio  is  still  increased  at  lower  turbine  inlet  temperatures  the 
exiting  steam  transitions  from  superheated  vapor  to  saturated  steam.  This  is  denoted  as 
pinch  point  location  #3  in  Figure  3.  With  further  increases  in  mass  ratio  the  pinch  point 
location  makes  one  last  transition  to  location  #4.  Throughout  the  temperature  range 
investigated  the  steam  always  exits  as  a  saturated  or  superheated  vapor.  At  higher 
turbine  inlet  temperatures  the  pinch  point  location  will  change  from  #6  to  #7  before 
moving  to  position  #4. 

At  higher  pressure  ratios  (16  and  20)  and  lower  turbine  inlet  temperatures  the 
maximum  efficiency  is  not  achieved  at  transition  from  pinch  point  location  #5  to  #6,  but 
at  the  transition  from  #3  to  #4.  Efficiency  increases  in  this  regime  when  the  pinch  point 
is  at  the  point  of  first  boiling.  In  this  area  increases  in  turbine  work  with  increasing  mass 
ratio  are  greater  than  the  penalties  of  heating  the  air  and  steam  to  the  turbine  inlet 
temperature.  The  penalty  of  heating  the  air  was  decreased  by  the  higher  compressor  exit 


temperatures  at  higher  compressor  pressure  ratios.  Figure  14  shows  plots  of  efficiency 
versus  steam/air  mass  ratio  for  increasing  turbine  inlet  temperatures  at  constant  pressure 
ratio  (Pr  =  16).  It  can  be  seen  here  that  the  maximum  efficiency  is  reached  at  higher 
mass  ratios  (the  transition  of  pinch  point  from  location  3  to  4)  at  the  lower  temperatures. 
As  the  temperature  increases,  less  and  less  gain  is  made  in  efficiency  as  mass  ratio 
increases.  The  plots  of  efficiency  versus  mass  ratio  get  very  flat  until  the  maximum 
efficiency  is  achieved  at  a  much  lower  mass  ratio  (the  transition  from  pinch  point  location 
#5  to  #6).  The  result  of  this  quick  transition  from  a  high  mass  ratio  to  a  lower  mass  ratio 
is  seen  in  Figures  12  and  13  where  the  steep  drop  in  specific  net  work  occurs. 

The  effect  of  limiting  the  steam/air  mass  ratio  to  preclude  any  condensation  in  the 
STIG  cycle  is  also  shown  in  Figures  4  through  13.  The  steam/air  mass  ratio  is  allowed 
to  increase  until  the  HRSG  stack  temperature  is  just  above  the  water  dewpoint 
temperature  for  the  steam's  partial  pressure.  In  general,  the  dewpoint  temperature 
increases  as  the  steam/air  mass  ratio  increases  or  as  the  steam  partial  pressure  increases. 
The  performance  curves  of  the  STIG  cycle  are  affected  when  a  high  steam/air  mass  ratio 
yields  the  maximum  efficiency.  This  occurs  at  the  higher  temperatures  for  the  lower 
pressure  ratio  curves  and  at  the  high  and  low  temperature  regions  of  the  higher  pressure 
ratios,  16  and  20.  The  STIG  performance  curves  with  the  dewpoint  temperature 
constraint  were  calculated  by  the  FORTRAN  program  called  CONDEN  found  in 
Appendix  B.  This  program  checks  the  STIG  performance  data  files  calculated  by  the 
main  program  GT.  If  the  stack  temperature  at  the  maximum  efficiency  is  below  the 
dewpoint  temperature,  the  steam/air  mass  ratio  is  decreased  until  the  a  stack  temperature 
is  achieved  which  is  above  the  dewpoint  temperature.  Program  CONDEN  uses  a 
relation  for  steam  saturation  temperature  as  a  function  of  saturation  pressure  given  by 
Irvine  and  Liley  [23]. 

Selected  performance  data  calculated  by  this  study's  program  are  compared  to 


data  at  the  same  parameters  presented  in  Bhutani  et  al.  [9]  and  Boyle  [11].  Parameters  of 
temperature  inlet  temperature,  steam/air  mass  ratio,  compression  pressure  ratio, 
isentropic  component  efficiencies,  ambient  conditions,  and  HRSG  pinch  point 
temperature  difference,  were  matched  to  those  of  the  literature  and  the  performance  char¬ 
acteristics  were  computed.  Table  4  shows  comparison  of  data  with  that  of  Bhutani  et  al. 
[9].  Table  5  shows  comparison  of  data  with  that  of  Boyle  [11],  As  can  be  seen  there  is 
close  agreement  between  the  study's  efficiency  data  and  that  of  the  literature.  There  is 
more  difference  in  the  results  of  net  specific  work.  These  are  due  to  different  methods  in 
pressure  loss  modeling  and  how  the  working  fluid  specific  heat  (Cp)  is  found. 

Second  Law  Analysis 

By  conducting  a  Second  Law  analysis  of  the  performance  of  the  Brayton  cycles 
further  insight  was  sought  into  the  performance  characteristics  of  these  cycles.  Second 
Law  efficiency  expresses  heat  inputs  as  Carnot  equivalent  work  requirements.  This 
means  that  for  a  certain  amount  of  heat,  one  could,  at  best,  only  utilize  the  Carnot 
equivalent  of  that  heat  input.  This  assesses  smaller  heat  input  penalties  in  calculating  the 
Second  Law  efficiency.  A  Second  Law  analysis  also  accounts  for  the  availability  of  inlet 
and  exit  streams.  This  can  assess  further  penalty.  In  general.  Second  Law  efficiencies 
are  higher  than  First  Law  efficiencies.  But  the  Second  Law  efficiencies  do  not  mirror 
First  Law  efficiency  in  most  of  the  cases  studied  here. 

Figure  15  shows  Second  Law  efficiency  versus  turbine  inlet  temperature  plots  for 
the  GT  simple  Brayton  cycle  at  the  three  highest  pressure  ratios  (12,  16,  20).  Figure  16 
plots  the  relation  of  the  same  parameters  for  the  GT  cycle  but  this  time  the  Second  Law 
efficiency  is  calculated  without  the  availability  of  the  exit  stream.  As  was  mentioned  in 
Chapter  II,  accounting  for  the  availability  of  this  stream  may  not  have  useful  purpose  as 
it  is  exhausted  out  the  stack  and  is  not  really  available.  For  the  case  of  the  GT  cycle  there 
is  a  noticeable  difference  in  whether  it  is  accounted  for  in  the  denominator  or  not.  When 
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the  exhaust  stream  availability  is  deducted  from  the  penalty  side  of  the  Second  Law 
efficiency,  the  efficiency  continues  to  increase  with  turbine  inlet  temperature.  When  the 
exhaust  stream  availability  is  not  accounted  for  there  is  a  maximum  Second  Law 
efficiency  within  the  studied  temperature  range.  This  is  due  to  the  effect  of  converting 
the  heat  input  to  a  Carnot  equivalent  work  term  which  asymptotically  approaches  a  limit 
as  temperature  increases.  This  provides  the  insight  that  when  expressing  heat  penalties 
as  Carnot  equivalent  energy  values,  the  GT  cycle  has  a  temperature  limit  for  maximum 
efficiency. 

Figure  17  shows  the  same  set  of  plots  for  the  IGT  cycle.  Here  the  dashed  lines 
are  the  Second  Law  efficiency  plots  calculated  without  accounting  for  the  exhaust  stream 
availability.  There  is  not  much  difference  between  these  plots,  whether  this  availability 
is  considered  or  not.  This  is  due  to  the  small  quantity  this  availability  represents. 
Because  this  system  has  a  regenerator  which  takes  advantage  of  the  hot  turbine  exhaust 
to  increase  efficiency,  the  cycle's  final  exhaust  temperature  is  much  lower  than  the  GT 
cycle  at  comparable  conditions.  The  exhaust  stream  availability  is  then  much  lower  for 
the  IGT  cycle.  It  is  also  seen  that  there  is  no  maximum  Second  Law  efficiency  achieved 
for  the  IGT  cycle  over  the  studied  temperature  range.  The  effect  of  two  heating  stages 
does  not  provide  the  limiting  effect  of  Carnot  energy  representation  that  was  seen  for  one 
heating  stage.  For  this  cycle  the  Second  Law  efficiency  continues  to  increase  as 
temperature  increases. 

The  results  of  the  STIG  cycle  are  again  more  complex  than  the  other  two  systems. 
It  is  also  a  matter  of  the  performance  of  the  heat  recovery  steam  generator.  Three 
definitions  of  the  Second  Law  efficiency  for  steam  producing  cycles  were  calculated  for 
comparison.  The  first  efficiency  represents  the  STIG  cycle  as  depicted  in  Figure  1.  The 
second  two  efficiencies  represent  steam  produced  in  a  HRSG  without  injection  into  the 
gas  turbine.  These  two  cases  represent  steam  produced  for  cogeneration  purposes.  On 


the  first  of  these  two  cases,  the  steam  is  heated  to  the  turbine  inlet  temperature.  In  the 
second  case  the  steam  is  not  heated  after  it  leaves  the  HRSG.  To  handle  these  last  two 
cases,  the  study’s  program  must  calculate  the  turbine  and  HRSG  processes  with  only  air 
as  the  working  fluid. 

To  more  easily  represent  what  the  Second  Law  analysis  predicts  about  the  STIG 
cycle  performance,  efficiency  plots  are  calculated  versus  steam/air  mass  ratio  at  a  number 
of  turbine  inlet  temperatures.  This  shows  the  effect  of  mass  ratio  and  turbine  inlet 
temperature  on  performance.  The  effect  of  pressure  ratio  is  not  as  significient  as  these 
other  two  parameters  but  will  be  discussed.  First  Figure  18  shows  the  effect  of  mass 
ratio  on  First  Law  efficiency  for  a  range  of  temperatures  at  a  pressure  ratio  of  16.  It  has 
been  explained  in  detail  earlier  how  the  maximum  efficiency  is  achieved. 

For  the  plots  of  Second  Law  efficiency  for  the  STIG  cycle  it  was  determined  to 
calculate  the  efficiency  without  accounting  for  the  availability  of  the  exit  stream.  Figure 
19  shows  the  difference  between  the  two  calculations.  The  dashed  line  is  the  efficiency 
without  the  exit  stream  availability  accounted  for.  At  low  mass  ratios  it  is  such  a  large 
term  compared  to  the  heat  penalties  that  it  shows  that  it  is  better  not  to  use  any  steam. 

Figure  20  shows  the  effect  of  mass  ratio  on  the  first  definition  of  Second  Law 
efficiency.  This  represents  the  Second  Law  efficiency  for  the  STIG  system  depicted  in 
Figure  1.  It  can  be  seen  that  a  maximum  efficiency  is  not  achieved  at  lower 
temperatures.  A  maximum  efficiency  is  later  achieved  as  the  turbine  inlet  temperature  is 
increased.  At  lower  turbine  inlet  temperatures  the  Second  Law  efficiency  continues  to 
increase  when  the  pinch  point  is  at  locations  6,  3,  or  4.  The  difference  in  the  Second 
Law  efficiency  and  the  First  Law  efficiency  in  this  regime  is  that  for  the  Second  Law 
analysis  the  heat  penalty  side  is  not  increasing  as  fast  as  the  net  work  side  of  the 
efficiency  calculation.  For  the  First  Law  calculation  the  opposite  is  true  in  this  regime. 
The  reason  for  the  dropping  off  of  the  heat  penalty  for  the  Second  Law  efficiency  can  be 
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Figure  19.  Second  Law  efficiency  versus  steam/air  mass  ratio  for  STIG  cycle  at 
turbine  inlet  temperature  of  1500  K  fPr  =  16). 
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found  in  its  method  of  calculation.  The  effect  of  increasing  turbine  inlet  temperature  is  to 
decrease  the  heat  penalty  for  heating  the  steam  to  superheat,  when  the  steam  exits  the 
HRSG  as  saturated  vapor,  with  pinch  point  locations  at  #3  and  #4.  Increasing  turbine 
inlet  temperature  also  decreases  the  heat  penalty  to  heat  the  steam  to  the  desired  turbine 
inlet  temperature,  when  the  steam  exits  the  HRSG  as  superheated  vapor,  with  pinch 
point  location  at  #6.  This  results  in  a  wider  temperature  range  than  the  First  Law 
analysis,  where  the  efficiency  increased  for  pinch  point  locations  at  points  6  and  3. 
There  is  a  limit  to  the  effect  of  decreasing  these  heat  penalties  by  increasing  turbine  inlet 
temperature.  With  increasing  turbine  inlet  temperature,  the  transition  to  pinch  point 
locations  6  and  3  happens  at  higher  mass  ratios.  With  a  higher  mass  ratio,  the  effect 
previously  mentioned  is  reversed  by  the  higher  quantity  of  steam. 

Increases  in  pressure  ratio  have  a  similar  effect  in  increasing  the  temperature  range 
where  Second  Law  efficiency  continues  to  increase  as  steam/air  mass  ratio.  For 
increasing  pressure  ratio,  the  heat  penalty  to  heat  the  exiting  steam  through  superheat  to 
the  desired  turbine  inlet  temperature  decreases,  plus  the  heating  penalty  of  heating 
superheat  steam  decreases.  In  the  case  of  pinch  point  #6,  the  exiting  superheat  steam 
increases  in  temperature  as  the  pressure  ratio  increases.  For  pinch  point  locations  #3  and 
#4,  the  saturation  temperature  at  which  the  steam  exits  increases  as  the  pressure  ratio 
increases.  These  increases  in  exit  temperature  decrease  the  required  steam  heat  inputs  as 
mentioned  above. 

In  total,  the  analysis  of  the  Second  Law  efficiency  shows  that  at  lower  turbine  inlet 
temperatures  the  efficiency  increases  throughout  the  steam/air  mass  ratio  range.  At 
progressively  higher  temperatures,  as  pressure  ratio  increases,  a  maximum  Second  Law 
efficiency  is  finally  achieved  at  the  same  steam/air  mass  ratio  as  the  maximum  First  Law 
efficiency. 

For  comparison  purposes,  a  Second  Law  analysis  is  done  on  a  GT  cycle  which 


produces  steam  for  cogeneration  purposes  but  has  no  steam  injection.  Figure  21  shows 
similar  Second  Law  efficiency  versus  steam/air  mass  ratio  plots  at  various  turbine  inlet 
temperatures  for  the  system  which  heats  the  steam  to  the  turbine  inlet  temperatures.  This 
figure  shows  that  the  Second  Law  efficiency  continues  to  increase  throughout  the  mass 
latio  range  for  all  the  selected  temperature  plots.  A  maximum  efficiency  is  reached  at  a 
higher  turbine  inlet  temperature  (not  shown  in  Figure  21)  than  any  of  the  other  systems 
evaluated. 

This  system  experiences  the  same  decreases  in  some  of  the  heat  penalties  at  pinch 
point  locations  #6  and  #3  as  explained  for  the  previous  systems.  In  addition,  this  system 
accounts  for  the  availability  of  the  steam  at  the  turbine  inlet  temperature  on  the  net  work 
side  of  the  efficiency  ratio.  This  has  a  positive  effect  on  the  Second  Law  efficiency  as 
turbine  inlet  temperature  increases.  By  accounting  for  the  steam  in  this  manner,  the 
Second  Law  efficiency  continues  to  increase  with  mass  ratio  over  a  wider  temperature 
range  than  the  STIG  cycle  and  the  system  to  be  shown  next. 

Figure  22  shows  the  Second  Law  efficiency  versus  steam/air  mass  ratio  plots  for 
various  temperatures  at  Pr  =  16  for  the  system  in  which  the  steam  is  not  heated  after 
exiting  the  HRSG.  For  this  system  there  is  always  a  maximum  Second  Law  efficiency 
achieved  within  the  mass  ratio  and  temperature  ranges  of  this  study.  At  lower  turbine 
inlet  temperatures,  the  maximum  efficiency  occurs  at  the  transition  from  pinch  point  #3 
to  #4.  At  this  point  the  availability  of  the  steam  exiting  the  HRSG  decreases  with  the 
pinch  point  located  at  position  #4.  At  higher  turbine  inlet  temperatures  the  maximum 
efficiency  occurs  at  lower  mass  ratios  where  the  pinch  point  location  transitions  from  #5 
to  #6.  The  system  previous  to  this  showed  that  higher  efficiencies  can  be  achieved  by 
further  heating  the  steam  from  the  point  left  off  at  by  this  system. 

As  a  note,  in  calculating  the  performance  of  a  HRSG  with  only  air  in  the  hot 
stream,  the  exiting  steam  can  leave  the  HRSG  as  liquid  water,  saturated  steam,  or 


temperatures  (Pr  =  16). 


0.75 


HRSG  where  steam  is  not  heated  beyond  the  HRSG  at  five  turbine 
inlet  temperatures  (Pr  -  16). 


CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  First  Law  analysis  it  can  be  concluded  that  both  of  the  Brayton  derivative 
cycles,  IGT  and  STTG,  offer  significant  improvements  over  the  simple  GT  Brayton 
cycle.  In  this  study  a  major  emphasis  has  been  placed  on  examining  cycle  performance 
as  a  function  of  cycle  maximum  temperature.  No  attempts  are  made  to  define  equipment 
concepts  for  the  heat  addition  process.  As  turbine  inlet  temperature  is  increased,  all  three 
cycles  demonstrate  a  diminishing  increase  of  efficiency  as  temperature  increases  but  the 
IGT  and  STIG  are  quite  superior  to  the  GT  cycle.  The  IGT  cycle  is  the  clear  winner  in 
First  Law  efficiency  with  the  STIG  in  second.  The  STIG  cycle  shows  a  most  dramatic 
increase  in  net  work  output.  It  demonstrates  an  increasing  relation  with  temperature. 
The  STIG  work  output  can  be  more  easily  controlled  by  the  steam  injection  rate,  which 
need  not  require  fluctuations  in  rotating  machinery  speed.  The  STIG  cycle  represents  a 
much  simpler  system.  This  equates  to  smaller  capital  costs.  The  STIG  cycle  results 
shown  here  and  the  other  advantages  mentioned  in  Chapter  I  demonstrate  that  the  STIG 
cycle  is  the  most  promising  candidate  for  higher  temperature  applications.  Table  6 
shows  the  improvement  in  the  net  work  output  and  efficiency  that  the  IGT  and  STIG 
cycles  have  over  the  GT  cycle  at  selected  maximum  temperatures  and  a  constant  pressure 
ratio. 

As  pointed  out  in  Chapter  II  the  limits  on  turbine  inlet  temperature  due,  to  the 
effectiveness  of  turbine  cooling  schemes  and  the  turbine  material's  temperature  limit,  do 
not  allow  application  of  the  whole  range  of  the  temperatures  investigated  here.  Further 
technological  advances  will  have  to  be  made  before  the  application  of  the  higher 


68 


temperatures  studied  here  can  be  realized.  Within  current  temperature  limits  the  STIG 
cycle  still  exhibits  a  greater  net  work  capacity  although  its  greatest  improvements  over 
the  IGT  cycle  in  net  work  are  achieved  at  higher  temperatures.  The  STIG  is  not  as 
efficient  as  the  IGT  in  this  range  but  offers  several  other  improvements  of  simplicity, 
lower  capital  costs,  and  flexibility  in  the  control  of  the  steam  injection  to  vary  work 
output  or  turbine  inlet  temperature.  These  advantages  are  deemed  more  valuable  overall 
than  the  IGT's  improvement  in  efficiency.  And  with  further  advances  in  maximum 
turbine  inlet  temperature,  the  STIG  cycle  will  offer  greater  improvements  in  net  work. 

The  Second  Law  analysis  provided  some  further  insights  in  the  performance  of 
these  cycles.  When  the  efficiency  is  represented  as  the  ratio  of  net  work  output  to 
reversible  work  the  GT  cycle  experienced  maximum  efficiencies  at  intermediate  turbine 
inlet  temperatures.  The  Second  Law  efficiencies  calculated  for  the  IGT  cycle  was  similar 
to  its  First  Law  efficiencies.  Its  Second  Law  efficiency  continued  to  increase  with 
turbine  inlet  temperature.  Like  the  First  Law  efficiency,  the  increase  in  the  Second  Law 
efficiency  diminished  in  the  higher  temperature  range.  The  maximum  Second  Law 
efficiency  for  the  STIG  cycle  was  achieved  at  a  higher  steam/air  mass  ratio  than  the  First 
Law  efficiency  at  the  lower  temperature  ranges  for  each  pressure  ratio.  At  the  higher 
pressure  ratios  no  maximum  efficiency  was  achieved  within  the  steam/air  mass  ratio 
range  at  these  lower  temperatures.  This  demonstrates  that  for  the  STIG  cycle,  less 
irreversibilities  are  experienced  at  higher  mass  ratios  than  those  which  achieved  the 
maximum  First  Law  efficiencies  at  these  lower  temperatures.  Location  of  the  higher 
irreversibilities  within  the  cycle  was  not  included  in  the  scope  of  this  study. 

Both  First  and  Second  Law  efficiencies  have  been  provided  for  each  cycle  and  for 
a  GT  cycle  with  a  HRSG.  Comparisons  are  made  of  when  maximum  efficiencies  are 
achieved  as  the  parameters  of  turbine  inlet  temperature,  steam/air  mass  ratio,  and 
compressor  pressure  ratio  are  varied.  But  for  solar  power  generation  applications  which 
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is  more  applicable,  First  or  Second  Law  efficiency?  The  utility  of  the  long  used  First 
Law  efficiency  is  its  value  of  measuring  what  you  are  getting  for  your  money.  The  First 
Law  efficiency  tells  the  plant  manager  what  amount  of  power  is  generated  for  the 
associated  fuel  costs.  For  a  solar  powered  system  there  are  no  fuel  costs.  The  Second 
Law  efficiency  makes  more  sense  in  telling  the  designer  how  close  a  given  design  is  to  a 
completely  reversible  system.  And  with  a  more  detailed  Second  Law  analysis  of  each 
system  component  the  designer  can  locate  where  the  major  irreversibilities  are. 

Where  the  parameters  of  this  study  provide  a  maximum  in  First  Law  efficiency  but 
not  in  Second  Law  efficiency,  it  would  be  more  useful  to  solar  power  applications  to 
seek  the  parameters  which  provide  the  maximum  Second  Law  efficiency  as  shown  in 
this  study. 

Each  of  the  simplifying  assumptions  of  this  study  represents  recommended  areas 
of  further  study  in  this  area.  There  are  two  recommendations  which  should  be  among  a 
more  detailed  study  of  these  power  cycles.  First,  a  turbine  cooling  scheme  should  be 
used  that  varies  the  amount  of  compressor  air  bleed  off  with  the  turbine  inlet 
temperature.  Bhutani  et  al.  [9]  made  a  good  case  for  the  straight  5%  bleed  off  of 
compressor  air  at  a  turbine  inlet  temperature  of  1700  K.  This  study  used  this  method. 
To  provide  a  conservative  performance  analysis  for  higher  temperatures,  a  cooling 
system  which  tasks  more  air  for  turbine  cooling  may  be  better. 

A  second  recommendation  for  further  study  is  to  adopt  a  set  pressure  loss  figure 
for  each  applicable  process  from  among  the  literature.  Table  7  lists  the  pressure  losses 
calculated  for  the  STIG  and  IGT  cycles  at  various  parameters.  The  method  developed  in 
this  study  for  calculating  the  pressure  loss  in  terms  of  average  temperature  and  initial 
pressure  is  applicable  for  only  selective  ranges  of  these  two  parameters.  This  judgement 
is  made  by  comparing  the  results  in  Table  7  to  selected  pressure  loss  figures  from  the 
literature  listed  in  Table  8. 


iaoic  / 

Pressure  Losses  Calculated  in  This  Studya 


Cycle 

Component 

Pr 

%Apb 

STIGC 

Heater 

12 

0.9 

16 

0.5 

HRSG 

12 

0.7 

16 

0.4 

IGT 

Intercooler 

12 

3.6 

16 

2.8 

Regenerator 

12 

0.7 

16 

0.4 

First  Heating 

12 

1.2 

16 

0.6 

Second  Heating 

12 

13.9 

16 

10.1 

^Turbine  inlet  temperature  for  each  example  is  1500  K. 

b%AP  represents  percent  pressure  drop  (%AP  =  APUOOyPjnjtjai). 

cMr  for  STIG  cycle  examples  is  0.250. 


Table  8 

Pressure  Loss  Schemes  Listed  in  Literature 


Reference 

%APa 

Comments 

Fraize  and  Kinney  [5] 

2.0 

Assessed  at  each  heating  process  and  HEX. 

Bhutani  et  al.  [9] 

4.0 

Assessed  at  combustor  only. 

Boyle  [ll]b 

Assessed  at  burner 

Assesed  to  turbine  exhaust  flow  in  HRSG 

Assessed  to  steam  flow  in  HRSG 


There  are  numerous  further  adaptations  and  more  detailed  modeling  schemes 
which  are  applicable  to  more  advanced  study  in  this  area.  Studying  the  effect  of  changes 
in  ambient  conditions  would  not  require  any  adaptation  to  the  existing  program.  A 
further  change  to  the  STIG  cycle  without  an  addition  of  rotating  machinery  is 
intercooling  to  reduce  compressor  work.  Where  this  study  left  off  in  its  Second  Law 
analysis,  further  studies  can  analyze  individual  component  irreversibilities.  There  is 
certainly  more  to  be  done  in  this  area.  A  goal  of  this  study  is  to  present  its  findings  in  a 
manner  to  give  further  work  in  this  area  several  possible  stepping  off  points. 


PROGRAM  GT 


Program  finds  1st  and  2d  law  cycle  efficiency  and  net  work  versus 
turbine  inlet  temp  for  simple  Brayton  cycle  (GT);  Brayton  cycle 
with  intercooling  ,  reheat,  and  regeneration  (IGT);  and  Brayton 
cycle  with  steam  injection  (STIG). 

CHARACTERS  CYCLE,  OUTPUT 
INTEGER  C,I,J,K,K1 ,L,LAW,II,JJ,R 

REAL  T6,T(9),P(9),Te(9),Pe(9),T5s,T5se,ErA,EPS(6),Wc(2),Wct, 

+Wt ( 2 ) , Wtt , Wp , Wpe , Wee , Wte , Pro ( 5 ) , PL , Mr , Cpmix , Cph , Cpsc , Qha( 2 ) , Qhs , 
+Qh , Qhex , Qhat , Qhe , Qhse , Qhexe , Qhat  e , Cpmixe , Cphe , Cpsce , Ke , Tstack , 
+Tstace,MAX(1334,5),MAX1 (491 ,4) ,Wnet,Wnete,Be,We(3) ,Wt1 ,Bs1  ,Bs2, 

+Be1 ,Wes(3),T5s1 ,Tstak1 ,Qhexa,T9a 

DATA  Pro ( 1 ) /4 • 0/ , Ke/ . 4299 / , SL41  /4 • 28571 42E-04/ , SL42/4 . 31 03448E-04/ , 
+SL81  /4- 1 1 21 495&-04/,  SL82/3. 960396E-04/,  SL1 21  /3-  8709677E-04/, 

+SL1 22/3 . 771 9298E-04/,  SL1 61/3-5251 798E-04/ ,  SL1 62/3 . 6885245E-04/, 
+SL201  /3.8E-04/,SL202/3*8636363E-04/ 

Make  choice  of  which  cycle  to  be  analyzed. 

¥RITE(*,10) 

FGRMAT( '0' , 'Enter  code  below  for  desired  cycle  analysis:'/'  ',T5, 

+ ' G-For  simple  Brayton  cycle' /'  ' , T5 , ' I-Eor  improved  Brayton  cycle' 
+, '  with  intercooling, regene rat ion,&reheat ' /'  ',T5,'S-Por  Brayton  ' 
+,' cycle  with  steam  injection') 

READ(*f 20)  CYCLE 
FORMAT (A1 ) 

Make  choice  of  which  compression  pressure  ratio  to  be  analyzed. 
WRITE(*,30) 

FORMAT ( 'O' . 'Enter  integer  below  for  desired  compressor  pressure  ', 
+' ratio  (Pr):'/'  ’ ,T5, ' 1 -Pr=4 ' / '  ’  T5, ’2-Pr=8'/'  • ,T5, '3-Pr=12'/ 

+'  *,T5* * 4 — Pr=1 6 ' / '  ,,T5,'5-Pr=20') 

REAL (*,40)  L 
FQRMAT(II) 

Make  choice  of  which  type  of  output  is  desired. 

WRITE(*,50) 

FGRMAT( 'O' , 'Enter  responce  for  desired  output  below:'/'  ',T5, 
+'T-Tabular  data  output'/'  ',T5,'1-Plot  of  turbine  inlet  ', 
+'temperature  vs  efficiency'/'  ',T5,'2-Plot  of  steam/air  mass  ', 

+' ratio  vs  efficiency'/'  ',T5, '5-Plot  of  First  Law  efficiency  ', 
+'and  corresponding  mass  ratio'/'  ' ,T7, 'versus  turbine  inlet  ', 

+' temperature' ) 

READ (*,60)  OUTPUT 
FORMAT (A1 ) 

Make  choice  of  which  type  of  efficiency  to  be  outputted. 


IF(OUTPUT.EQ. ’1 1 .OR. OUTPUT. EQ. '2' )  THEN 
WRITE(*,61 ) 

FCRMAT( 'O' , 'For  efficiency  plots  specify  which  type  of 
+ 'analysis  below: '/'  ' ,T5, '1 -First  Law  efficiency'/'  ',T5,’2-', 
+'Second  Law  efficiency') 

READ(*,62)  LAW 
F0RMAT(I1  ) 

MDIF 

Choose  which  turbine  inlet  temperature. 

IF(0UTPUT.NE.'1 ’.AND. OUTPUT. NE. '3')  THEN 
WRITE(*,63) 

F0RMAT( 'O' , 'Enter  value  of  turbine  inlet  temperature (K) ' , 

+  '(F6.1 ).') 

KEAD(*,65)  T6 
FORMAT (F6.1  ) 

ENDIF 

Choose  which  steam/air  mass  ratio  (Mr),  (for  STIG  cycle). 
IF(CYCLE.EQ. ’S')  THEN 

IF(OUTPUT.EQ. ' T '  .OR. OUTPUT. EQ.  '1  ' )  THEN 
WRITE(*,70) 

FORMAT ( 'O' , 'Enter  value  of  steam/air  mass  ratio  (Mr)', 

+'(F5.3).’) 

READ(*,80)  Mr 
FORMAT(F5-3) 

ENDIF 

ENDIF 

IF(CYCLE.NE.'S')  Mr=0.0 
IF(OUTPUT.NE.'T')  GO  TO  230 

For  tabular  data  output  for  specified  turbine  inlet  temperature 
and  specified  Mr. 

CALL  CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp,Wnet,Cph,Cpsc, 
-+Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1 ,Bia,Biw,Be1  ,Wes,K1 , 
+T5s1 ,Tstak1 ,Qhexa,T9a»Bs1  ,Bs2) 

DO  90  1=1 ,9 

Te(I)=1 .8*T(I) 

Pe(I)=0.14504*P(I) 

CONTINUE 
Wct=Wc(1 )+Wc(2) 

Wtt=Wt(1  )+Wt(2) 

Qhat=Qha(1 )+Qha(2) 

Qh=Qha(1  )+Qha(2)+Qhs 

Wce=Ke*Wct 

Wpe=Ke*Wp 

Wte=Ke*Wtt 

Wnete=Ke*Wnet 

Qhate=Ke*Qhat 

Qhse=Ke*Qhs 


Qhe=Ke*Qh 
Qhexe=Ke*Qhex 
Cpsce= . 2388*Cpsc 
Cpmixe=. 2388*Cpmix 
Cphe=.2388*Cpb 
T5se=1 .8*T5s 
Tstace=1  .8*Tstaek 

IF(CYCLE.NE. 'I')  THEN 
IF(CYCLE.EQ.'G')  THEN 
WRITE (1 5, 100)  PL 

FORMAT  ('O',  'SIMPLE  BRAYTON  CYCLE  WITH  OVERALL  Pr=  \F4.1) 
GO  TO  120 
ENDIF 

WRITE(1 5»  1 10)  PL, Mr 

FORMAT ('O', 'BRAYTON  CYCLE  WITH  STEAM  INJECTION,  Pr=  \F4.1, 
+2X, 'Mr=  ' ,F5-3) 

WRITE(15,130) 

FORMAT ('O', 'STAGE',  4X, 'T(K) ' ,6X, ’T(R) ’ ,5X, 'P(kPa) ' ,3X, 
+'P(psi) '/'O' ) 

1=1 

WRITE(15,140)  I,T(1  ),Te(1  ),P(1  ),Pe(1 ) 

F0RMAT('  ' ,2X,I1 ,4X,F8.3,2X,F8.3,2X,F8.3,2X,F8.3) 

1=2 

WRITE(1 5,140)  I,T(5),Te(5),P(5),Pe(5) 

IF(CYCLE.EQ. 'S')  THEN 
WRITE(1 5, 1 50)  T5s,T5se 

FORMAT ('  ', 'STEAM  IS  INJECTED  INTO  COMBUSTOR  AT  \F6.1, 
+'K(',F6.1,'R)') 

ENDIF 

1=3 

WRITE (15, 140)  I,T(6),Te(6),P(6),Pe(6) 

1=4 

WRITE(15,140)  I,T(9),Te(9),P(9),Pe(9) 

IF(CYCLE.EQ.'S')  THEN 

WRITE(15,160)  Tstack,Tstace,K,K1  ,T5s1  ,Tstak1  ,Qhexa,T9a 
FORMAT ( '  ' , 'EXHAUST  STACK  TEMPERATURE  IS  ' ,F6. 1 , 'K( ' , 
+F6.1,'R)'/'  ','HRSGmode  is  ',11/'  ', 'Air  HRSG  Results  are: ' , 

+'  Mode  is  ',11/'  ',22X, 'Water  leaves  at  temp  ' , F6 - 1 / '  ',22X, 
+'Tstack  is  ' ,F6. 1 /’  \22X,'Qhex  is  ',F8.3/'  ' ,22X, 'Exhaust  ', 

+' enters  HRSG  at  ',F6.1) 

BIDIF 

WRITE( 15,170) 

FORMAT( 'O'/’O' , 10X, ' (kJ/kg) ’ , IX, ' (BTU/lbm) ' ) 

WRITE ( 1 5 , 1 80 )  Wet , Wee , Wp , Wpe , Wtt , Wte , Wnet , Wnete , Qhat , Qhate , 
+Qhs,Qhse,Qh,Qhe,Qhex,Qhexe,Be,We(1 ),We(2),We(3) 

FORMAT ( '  ', 'Wc',8X,F7.3,2X,F7.3/'  ', 'Wp',8X,F7.3,2X,F7.3/' 
+'Wt',7X,F8.3,2X,F7.3/'  ' , 'Wnet' ,5X,F8.3,2X,F7.3/’  ','QHair', 
+4X,F8.3,2X,F7.3/'  ' , 'QHsteam' ,2X,F8.3,2X,F7.3/'  ','QHtotal', 
+2X,F8.3,1X,F8.3/’  ’ , 'Qhrsg' ,4X,F8.3,2X,F7.3/’  ' , 'Bexhaust' , 
+1X,F8.3/'  ','We(1 )',4X,F8.3/'  ' , ’We(2) ' ,4X,F8.3/'  ' , ’We(3) ' ,4X, 
+F8.3/'0') 

IF(CYCLE.EQ. 'S')  THEN 

WRITE(1 5,185)  Wtl ,Bia,Biw,Be1 ,Wes(1 ),  'es(2),Wes(3),Bs1 ,Bs2 
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For  data  output  of  cycle  efficiency  and  specific  net  work 
versus  steam/air  mass  ratio  (0.01-0.50),  (3TIG  cycle). 

IF (OUTPUT. EQ. '2' )  THEN 
Mr =0.01 

DO  270  JJ=1 ,491 

CALL  CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp,Wnet,Cph, 
+Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1 ,Bia,Biw,Be1 , 
+Wes,K1  ,T5s1  ,Tstak1  ,Qhexa,T9a,Bs1  ,Bs2) 

IF(LAW.EQ.I)  WRITE(1 7,255)  Mr , ETA , Wnet , K 
FORMAT ( '  ',F5.5,2X,F6.4,2X,F8.3,2X,I1 ) 

IF(LAW.EQ. 2)  WRITE(17,260)  Mr,(EP3(I),  I=1,6),K,K1 
FORMAT ('  ',F5*3,2X,6(F6.4,2X),I1,2X,I1  ) 

Mr=Mr+0.001 

CONTINUE 

ENDIF 
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For  output  of  maximum  First  Law  efficiency  and  corresponding  mass 
ratio  &  specific  net  work  versus  turbine  inlet  temperature  (STIC 
cycle ) . 

IF( OUTPUT. EQ. '3 ' )  THEN 
T6=1 000.0 
IF(L.EQ.I)  R=823 
IF(L.EQ.2)  R=1 050 
IF(L.EQ.3)  R=1 1 78 
IF(L.EQ.4)  R=1 266 
IF(L.EQ.5)  R=1 334 
DO  310  11=1  ,R 

IF ( T6 . LT . 1 073 • • AND . L. EQ . 4 )  THEN 
Mr=0.01 0 
C=491 
GO  TO  275 
ENDIF 

-IF(T6.LT.  1 1 60. . AND. L. EG. 5)  THEN 
IF(T6.LT. 1 066. )  THEN 
Mr=0.414 

CALL  CYCLE2 ( T6 , L , Mr , CY OLE , T , P, PL, ETA , EPS , Wc , Wt , Wp . 
+Wnet , Cph , Cpsc , Cpmix , Qha, Qhs , Qhex , Ts tack , T5s , K, Be , We , Wtl  , Bia, Biw , 
+Be1 ,Wes,K1 ,T5s1 ,Tstak1  ,Qhexa,T9a,Bs1  ,Bs2) 

MAX (1 1, 1  )=T6 
MAX(II,2)=Mr 
MAX(II,3)=ETA 
MAX(II,4)=Wnet 
MAX (11,5) =Tstack 
GO  TO  295 
ENDIF 
Mr=0.01 0 
C=491 
GO  TO  275 
ENDIF 

IP(L.EQ.I)  THEN 

Mr=SL41  *(T6-1 000.  )-K).  1 3 
C =(  3142  *  ( T 6- 1 000 .)  -hj  .  1 6-M  r )/.  001  + 1  . 0 


m.  i 
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IF(L.EQ.2)  THEN 

Mr=SL81  *(T6-1000. )+0.06 
C=(SL82*(T6-1000.  )+0.10-Mr)/. 001+1 .0 
ENDIF 

IF(L.EQ.3)  THEN 

Mr=SL121*(T6-1000. )+0.02 
C=(SL122*(T6-1000. )+0.05-Mr)/. 001+1 .0 
ENDIF 

IF(L.EQ.4)  THEN 

Mr=SL1 61 *(T6-1 000. )+0.01 
C=(SL162*(T6-1000.)40.05-Mr)/. 001+1 .0 
ENDIF 

IF(L.EQ.5)  THEN 

Mr=SL201 *(T6-1 1 50. )+0.025 
C=(SL202*(T6-1 1 50. )+0.075-Mr)/. 001+1 .0 
ENDIF 

275  DO  280  JJ=1 ,C 

CALL  CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,¥c,¥t,¥p,¥net, 
+Cph,Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,¥e,¥t1 ,Bia,Biw,Be1 , 
+¥es,K1  ,T5s1  ,Tstak1 ,Qhexa,T9a,Bs1 ,Bs2) 

MAXI (JJ,1  )=Mr 
MAXI (JJ,2)=ETA 
MAXI (JJ,3)=Wnet 
MAXI  (JJ,4)=Tstack 
Mr=Mr+0.001 
280  CONTINUE 

MAX(II,3)=MAX1 (1,2) 

DO  290  1=2  C 

IF(MAX1(I,2).GT.MAX(II,3))  THEN 
MAX(II,1 )=T6 
MAX(II,2)=MAX1 (1,1 ) 

MAX(II,3)=MAX1 (1,2) 

MAX(II,4)=MAX1 (1,3) 

MAX(II,5)=MAX1 (1,4) 

ENDIF 

290  CONTINUE 

295  WRITE(18,300)  (MAX(II,J),  J=1,5),L 

300  FORMAT ('  \F6.1,2X,F5.3,2X,F6.4,2X,F8.3,2X,F9.5,2X,I1 ) 

T6=T6+1 .0 
310  CONTINUE 
INDIF 
C 

320  STOP 
END 
C 

C - SUBROUTINES - 

C 

c 

SUBROUTINE  CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,¥t,Wp,Wnet,Cph, 
+Cpec,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,¥e,¥t1 ,Bia,Biw,Be1 ,¥es, 
+K1  ,T5s1  ,Tstak1 ,Qhexa,T9a,Bs1  ,Bs2) 


C 

C 


This  subroutine  calculates  the  state  points  for  the  desired  cycle 


using  First  Law  energy  balances.  The  work  outputs/inputs  and  the 
heat  inputs  are  calculated.  Frcm  these  values,  the  First  and 
Second  Law  efficiencies  are  calculated. 

CHARACTER*  1  CYCLE,  CYCLE1 
INTEGER  I, J,K,K1  ,L 

REAL  T6,T(9),P(9),To5,Tn5,T5s,ETAc,ETAt,ETAp,EFA,Wc(2),Wt(2),Wp, 
+Kplost , Plost , Pr , Pro ( 5 ) , PL , Mr , Cpa , Cps , Cpmix , Cph , Cpsc , Qha( 2 ) , Qhs , 
+Qh,Qhex,RH0s,Tstack,Wnet,Cpa1 ,Cpa2,Cpa3,Cpa4,Cpm1 ,Tsat1 ,EPS(6), 
+Bia,Biw,Be,We(3),Ya,Ys,Te1  ,Te2,Cpm2,Cpm3,Cpm4,T9a,Wt1  ,Bs1  ,Bs2,Mr1 
4Ves(3),T5s1  ,Tstak1  ,Qsat 

DATA  Pro(1  )/4.0/,To5/0.0/,ETAc/0.er7/,ETAt/0.89/,ETAp/0.70/, 
+RH0s/997.0/,Kplost/1 2. 1 28/,MWa/28. 97/,MWs/l 8.01 5/,Ru/8. 31 434/, 
+Te/0.0/ 


T(1 )=300.0 

P(1  )=101 .325 

T(6)=T6 
DO  10  1=2,5 

Pro(I)=Pro(I-1 )+4-0 
CONTINUE 

P(9)=P(1 ) 

Pr=Pro(L) 

PL=Pro(L) 

IF(CYCLE.EQ.’I')  Pr=SQRT(Pr) 

Determine  required  compressor  specific  work  (Wc). 

CALL  CWQRK(T(1 ),P(1 ),Pr,ETAc,P(2),T(2),Wc(1 )) 

IF(CYCLE.NE. 'I' )  THEN 

T(5)=T(2) 

P(5)=P(2) 

ENDIF 

If  IGT  cycle  is  called  for,  find  results  of  intercooler 
and  2d  compressor. 

IF(CYCLE.EQ. ’I')  THEN 

Determine  temp  drop  and  pressure  lost  in  intercooler. 

T(3)=327.5 

Plost=( (T(1 )+T(2) )/(2.0*P(2) ) )*Kplost 
P(3)=P(2)-Plost 

Determine  required  work  in  2d  compressor. 

CALL  CWQRK(T(3),P(3),Pr,ETAc,P(4),T(4),Wc(2)) 

ENDIF 


Set  an  initial  guess  for  regenerator  outlet  temp,  T(5)* 


IF(CYCLE.EQ.  'I')  T(5)=T(4)+100.0 

If  analyzing  IGT  cycle,  the  program  returns  here  to 
iterate  the  effect  of  the  regenerator. 

Determine  pressure  lost  across  the  regenerator. 

IF(CYCLE.EQ. 'I')  THEN 

Plost=((T(4)+T(5))/(2.0*P(4)))*Kplost 

P(5)=P(4)-Plost 

ENDIF 

Determine  pressure  loss  across  combustor. 

Plost=((T(5)+T(6))/(2.0*P(5)) )*Kplost 
P(6)=P(5)-Plost 

IF(CYCLE.EQ.’I’)  P(7)=P(6)/SQHT(P(6)/P(9)) 

IF(CYCLE.NE. 'I')  P(7)=P(9) 

Determine  turbine  specific  work  (Wt). 

CALL  TWQRK(CYCLE,T(6),P(6),P(7),Mr,ETAt,T(7),Cpmix,Wt(1 )) 
IF(CYCLE.NE.'I')  T(9)=T(7) 

If  IGT  cycle  is  called  for,  find  results  of  reheat 
and  2d  turbine. 

IF(CYCLE.EQ. 'I')  THEN 
T(8)=T(6) 

Determine  pressure  loss  in  reheat. 

Plost=((T(7)+T(8))/(2.0*P(7)))*Kplost 
P(8)=P(7)-Plost 
Determine  2d  turbine  work. 

CALL  TWCRK(CYCLE,T(8),P(8),P(9),Mr,ETAt,T(9),Cpmix,Vt(2)) 
ENDIF 

Determine  steam  exit  temperature  from  HRSG,  required  heat  input 
to  raise  exiting  steam  to  turbine  inlet  temperature,  and  pump 
work. 

IF(CYCLE.EQ.  'S’)  THEN 

CALL  HRSG1  (T(1 ) ,T(6) ,T(9) ,Mr,L,Cph,Qhex,T5s,Cpsc,Qhs,Tstack,K 
■fTsat  1  f  Qsat ) 

Plost=((T(6)+T(1  ))/(2.0*P(5) ) )*Kplost 
Wp=Mr*(P(5)-P(1 )+Plost)/(RHOs*ETAp) 

EJDIF 

Determine  heat  input  to  air  in  combustor. 

CALL  CPAIR1 (T(5)tT(6),Cpa) 


Qha(1 )=0.95*Cpa*(T(6)-T(5)) 


82 


C 

C  If  IGT  cycle  is  called  for,  add  in  heat  input 

C  for  the  reheat. 

C 

IP( CYCLE. EQ. ' I  * )  THEN 

CALL  CPAIR1 (T(7),T(8),Cpa) 
Qha(2)=0.95*Cpa*(T(8)-T(7)) 

ENDIP 

C 

C  Iteration  to  account  for  regenerator  is  done  here. 

C 

IF(CYCLE.EQ. 'I')  THEN 
Tn5=T(5) 

IP(ABS(Tn5-To5).LT.0.1 )  GO  TO  30 
To5=Tn5 

T(5)=T(9)-50.0 
GO  TO  20 
ENDIP 

C 

C  Iteration  to  determine  Texhaust  from  regenerator. 

C 

30  IF(CYCLE.EQ. 'I')  THEN 

CALL  CPAIR1 (T(4) ,T(5) ,Cpa1 ) 

Te2=T(4)+100.0 

35  CALL  CPAIR1 (Te2,T(9),Cpa2) 

Te2=T(9)-Cpa1*(T(5)-T(4) )/Cpa2 
IF(ABS(Te2-Te1).GE.0.1)  THEN 
Tel =Te2 
GO  TO  35 
INDIF 
ENDIF 

C 

IF ( CYCLE. NE. 'I')  THEN 
Wc(2)=0.0 
Wt(2)=0.0 
Qha(2)=0.0 
DO  40  J=2,4 
P(J)=0.0 
T(J )=0.0 

40  CONTINUE 

DO  50  1=7,8 
P(I)=0.0 
T(I)=0.0 

50  CONTINUE 

ENDIF 

C 

IF(CYCLE.NE.'3')  THEN 
Wp=0.0 
Qhs =0.0 
Qhex=0.0 
Cph=0.0 
Cpec=0.0 
Tstack=0.0 


ENDIF 

Determine  First  law  cycle  efficiency  (ETA). 

Wnet=Wt(1 )+Wt(2)-Wc(1 )-¥c(2)-Wp 
ETA=Wnet/(Qha(1  )+Qha(2)+Qhs) 

Determine  Second  Law  efficiency  for  desired  cycle. 

IF(CYCEE.EQ. 'G' )  THEN 

CALL  CPAIR1 (T(1  ),T(9) ,Cpa1  ) 

CALL  CPAIR1 (T(5),T(6),Cpa2) 

Be=0.95*Cpa1 *(T(9)-T(1  )*(1  .O+L0G(T(9)/T(1 )))) 

We(1 )=0.95*Cpa2*(T(6)-T(5)-T(1  )*LOG(T(6)/T(5) ) ) 

EPS(1  )=Wnet/(We(1  )-Be) 

EPS(2)=Wnet/We(1 ) 

ENDIF 

IF(CYCLE.EQ. 'I' )  THEN 

CALL  CPAIR1  (T ( 1  ),T(9),Cpa1 ) 

CALL  CPAIR1 (T(5),T(6),Cpa2) 

CALL  CPAIR1  (T(7),T(8),Cpa3) 

CALL  CPAIR1 (T(3),T(8),Cpa4) 

Be=0.95*Cpa1  *(Te2-T(1  )*(1  .OLOG(Te2/T(1 ) ) ) ) 

We(1 )=0.95*Cpa2*(T(6)-T(5)-T(1  )*L0G(T(6)/T(5) ) ) 
We(2)=O.95*Cpa3*(T(8)-T(7)-T(1)*L0G(T(8)/T(7))) 
We(3)=O.95*Cpa4*(T(3)-T(2)-T(1)*L0G(T(3)/T(2))) 

EPS(1 )=Wnet/(We(1 )+We(2)+We(3)-Be) 

EPS(2)=Wnet/(We(1  )+We(2)+We(3) ) 

ENDIF 

IF(CYCLE.EQ. 'S' )  THEN 

CALL  CFM(T(1 ) .Tstack.Mr.Cpml  ) 

CALL  CPH20(T(1),T(6),Cps1) 

CALL  CPH20(T(1  ),T5s,Cps2) 

CALL  CPAIR1  (T(5),T(6),Cpa) 

Ysu=  (0. 95*MWs )  /  (0. 95*MWs+Mr*MWa) 

Ys=(Mr*M-rfa) /(0. 95*MWa-^r*MWa) 

Bia=T(1  )*Ru*LOG(  1 .0/Ya)/MWa 
Biw=Mr*T  ( 1  )*Ru*L0G(1 .0/Ys)/MWs 

Be=(Mr+0.95)*Cpm1*(Tstack-T(1  )*(1  .O+L0G(Tstack/T(1 )))) 

CYCLE1  ='G' 

Mr  1=0.0 

CALL  TWCRK(CYCLE1  ,T(6),P(6),P(9),Mr1  ,ETAt,T9a,Cpm2,Wt1 ) 
We(1)=0.95*Cpa*(T(6)-T(5)-T(1)*L0G(T(6)/T(5))) 
IF(K.EQ.3.0H.K.EQ.4)  THEN 
We(2 )  =  ( 1 . 0— T ( 1  )/Tsat1  )*Qsat 

We(3)=Mr*Cpec*(T(6)-T5s-T(1 )*L0G(T(6)/T5s) ) 

CALL  HRSG2(T(1 ),T(6),T9a,Mr,L,T5s1 ,Tstak1 ,Wes,K1 ,Qhexa,Bs1  ,Bs2) 
CALL  CPAIR1 (T ( 1 ),Tstak1 ,Cpa1 ) 

Bel =0.95*Cpa1 *(Tstak1 -T(1 )*(1 .0+LOG(Tstak1  /T(1 ) ) ) ) 

EPS(1 )=Wnet/(We(1 )4We(2)-+We(3)+Bia+Biw-Be) 


o  o  o  o  o  o 


EPS ( 2)=Wnet/(We ( 1  )+We ( 2 )+We ( 3 )+Bia+Biw) 

EPS(3)=(Wt1-Wc(1 )-Wp+Bs1 )/(We(1 )+Ves(1  )+Wes(2)4¥es(3)-Be1 ) 
EPS(4)=(Wt1-Wc(1  )-WpfBs1  )/(We(1  )+Wes(1  )+Wes(2)+¥es(3) ) 
EPS(5)=(Wt1-Wc(1  )-Wp+Bs2)/(We(1  )-Be1 ) 

EPS(6)=(Wt1-Wc(1 )-WpfBs2)/We(1 ) 

ENDIF 

RETURN 

END 


SUBROUTINE  CWQRK(T1 ,P1 ,Pr,ETAc, P2,T2,Wc) 

This  subroutine  finds  the  required  compressor  work  input. 
CHARACTERS  TYPE 

REAL  T1 ,P1 , Kao, Ka,Pr,P2,Ts2,ETAc,T2o,T2n,T2,Cpa,Ra,Tisen, Tact 

DATA  Kao/1 .397/, Ra/. 28700/ 

TYPE='COMP' 

Ka=Kao 
P2=P1  *Pr 

Ts2=Tisen(TYPE,T1 ,Pr,Ka) 

T2o=Tact(TYPE,T1 ,Ts2,ETAc) 

CALL  CPAIR1 (T1 ,T2o,Cpa) 

Ka=Cpa/ ( Cpa-Ra) 

Ts2=Tisen(TYPE,T1  ,Pr,Ka) 

T2n=Tact (TYPE, T1 ,Ts2,ETAc) 

IF(ABS(T2n-T2o) .GE.0.1 )  THEN 
T2o=T2n 
GO  TO  10 
ENDIF 
T2=T2n 

Wc=Cpa*(T2-T1 ) 

RETURN 

END 


SUBROUTINE  TW0RK(CYCLE,T1 , PI  ,P2,Mr,ETAt,T2,Cpmix,Wt) 

This  subroutine  finds  the  turbine  work  output. 

CHARACTER*!  CYCLE 
CHARACTER*4  TYPE 

REAL  T1 , PI ,P2,Mr,Kao,Ka,Kso,Ks,ETAt,Ra,Rs,T2,Cpmix,Wt,Ts2,Tisen, 
+T2o , T2n , Tact , Cpa , Cpao , Cpe , Cpso , Kmix , Pr 

DATA  Kao/1 -397/,Kso/l -327/,Cpao/1 .0035/,Cpso/1 . 8723 /,Ra/. 28700/, 


+Rs/.461 52/ 


TYPE='TURB' 

T2o=0.0 

Ka=Kao 

Ks=Kso 

Cpa=Cpao 

Cps=Cpso 

Pr=P1  /P2 

IF(CYCLE.NE.'S')  Kmix=Ka 

IP ( CYCLE . EQ. • S ' )  Kmix= (Mr*Cps+ . 95*Opa ) / ( ( Mr*Cps/Ks )+ ( . 95*Cpa/Ka ) ) 
Ts2=Tisen(TYPE,T1 ,Pr,Kmix) 

T2n=Tact(TYPE,T1 ,Ts2,ETAt) 

IF(ABS(T2n-T2o) .GE.0.1  )  THEN 
T2o=T2n 

CALL  CPAIR1 (T1 ,T2o,Cpa) 

CALL  CPH20(T1 ,T2o,Cps) 

Ka=Cpa/ ( Cpa-Ra ) 

Ks=Cps/ (Cps-Rs) 

GO  TO  10 
END  IF 
T2=T2n 

IP(CYCLE.NE. 'S' )  Cps=0.0 
Cjmix=(Mr*Cps+0.95*Cpa)/(0.95-hMr) 

Wt=Cpmix*( .95+Mr)*(T1-T2) 

RETURN 

END 


SUBROUTINE  HRSG1 (T1 ,T6,T9,Mr,L,Cph,Qhex,T5s,Cpec,Qhs,Tstack,K, 
+Tsat1  ,Qsat) 

This  subroutine  calculates  the  performance  of  the  HRSG  for  the 
STIG  cycle.  It  finds  the  exiting  steam  temperature,  the  exhaust 
gas  stack  temperature, and  the  heat  required  to  raise  the  exiting 
steam  to  the  turbine  inlet  temperature.  For  this  subroutine  the 
hot  exhaust  gas  is  a  mixture  of  steam  and  air. 

INTEGER  K,L 

REAL  Tsat(5),Hfg(5),Cpw(5),Q1 ,Q2,Q3,Qhex,Qhex1 ,Qhs,Tpp,Tpp1 ,Tpp2, 
+T1 ,T6,T9,Tstack,Tstak1 ,Tstak2,Tstak3,T5s,T5s1 ,T5s2,Cprnix,Cpmix1 , 
+Cpmix2,Cps,Cpsc,Cph,Ch1  ,Ch2,Mr,Tsat1 ,Qsat 

DATA  Tsat(1 )/417.3/,Tsat(2)/444.1/,Tsat(3)/461 .7/,Tsat(4)/475-2/, 
+Tsat(5)/486.2/,Hfg(1)/2132.4/,Hfg(2)/2046.2/,Hfg(3)/l984.0/, 
+Hfg(4)/1932.7/,Hfg(5)/l887.8/,Tpp1 /0.0/,T5s1/0.0/,Tstakl/0.0/, 
-KJpw(1)/4.205/,Cpw(2)/4.215/,Cpw(3)/4.225/,Cpw(4)/4.230/, 
•+Cpw(5)/4.245 / 

Tsatl =Tsat(L) 

Check  if  water  can  reach  saturation. 


Tpp=T 1+26.0 

Q1 =Mr*Cpw(L)*(Tsat(L)-T1 ) 

CALL  CPM(Tpp,T9,Mr,Cpmix) 

Cpmixl  =Cpmix 

Qhexl  = (Mr+0 . 95 ) *Cpmix* ( T9-Tpp ) 

CALL  CPH20(Tsat(L),T6,Cps) 

Cpsc=Cps 

IP (Qhexl .LT.Q1 )  THEN 

Qhs=Q1  -Qhexl  -+Mr*(Hfg(L)+Cpsc*(T6-Tsat  (L) ) ) 

Tstack=Tpp 

An  approximation  (due  to  Cpw  value)  of  water  temp  entering 
combustor. 

T5s=Qhex1 /(Mr*Cpw(L) )+T1 
IF((T9-T5s).LT.26.0)  THEN 

Pinch  point  is  at  approach  temperatures. 

T5s=T9-26.0 
Tstak2=T 1+50.0 
Qhex=Mr*Cpw(L)*(T5s-T1 ) 

IF(Qhex.GT.Q1 )  THEN 

T5s=Qhex1 /(Mr*Cpw(L) )+T1 
GO  TO  15 
ENDIP 

CALL  CFM(Tstak2,T9,Mr,Cpmix) 

Tstak2=T9-Qhex/( (Mr+O. 95 )*Cpmix) 

IF(ABS(Tstak2-Tstak1 ).GE.0.1 )  THEN 
Tstakl =Tstak2 
GO  TO  10 
ENDIF 

Tstack=Tstak2 

Cph=Cpmix 

Qhs=Mr*(Cpsc*(T6-Tsat(L))+Hfg(L)+Cpw(L)* 

+(Tsat(L)-T5s)) 

K=2 


The  pinch  point  location  is  at  the  approach  temperature 
location  and  the  water  exits  the  HRSG  as  a  liquid. 

GO  TO  90 
ENDIP 
Qhex=Qhex1 
Cph=Cpmix1 
K=1 


The  pinch  point  location  is  at  the  water  entry  point  and 
the  water  exits  the  HRSG  as  a  liquid. 

GO  TO  90 
ENDIP 


Check  if  steam  can  reach  superheat.  First  assume  pinch  point  is 


oooo  o  r\j  o  o  o  oooooo 
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C  at  entry  to  saturation  region. 

C 

17  Tpp=Tsat(L)+26.0 

Tstakl  =0.0 

20  Tstak2=Tpp-Q1  / ( (Mr-K). 95 )*Cpmix) 

IF(ABS(Tstak2-Tstak1  ).GE.0.1 )  THEN 
CALL  CPM(Tstak2,Tpp,Mr,Cpmix) 

Tstakl =Tstak2 
GO  TO  20 
ENDIF 

Tstak3=Tstak2 

CALL  CPM(Tstak2,T9,Mr,Cpmix2) 

Chi  =Cpmix1  *(Mr+0.95) 

Ch2=Cpmix2* ( Mr+0 . 95 ) 

Qhex=Ch2* (T9-Tstak2 ) 

Q2=Mr*(Hfg(L)+Cpw(L)*(Tsat(L)-T1 )) 

IF(Qhex.LT.Q2)  THEN 
Cph=Cpmix2 

Qhs=Q2-Qhex+Mr*Cpsc*(T6-Tsat(L) ) 

Tstack=Tstak2 
Qsat=Q2-Qhex 
K=3 

The  pinch  point  location  is  at  the  beginning  of  the  sat¬ 
uration  region  and  the  steam  exits  the  HRSG  in  saturation. 

Check  to  see  if  pinch  point  is  at  water  entry  point. 

IF( (Tstack-TI ).LT.26.0)  THEN 
Tstack=T 1+26.0 

Check  to  see  if  Mode  3  is  more  applicable. 

Tpp2=Tsat (L)+30. 0 
CALL  CPM(Tstack,Tpp2fMr,Cpmix) 

Tpp2=Q1  /( (Mr+0.95)'*'Cpmix)+Istack 
IF(ABS(Tpp2-Tpp1  ).GE.0.1 )  THEN 
Tppl  =Tpp2 
GO  TO  25 
ENDIF 

IF( (Tpp2-Tsat(L) ) .LT.26.0)  THEN 
GO  TO  17 
ENDIF 

Qhex=Ch1 *(T9-Tstack) 

Cph=Cpmix1 

Qhs=Q2-Qhex+Mr*Cpsc*(T6-Tsat(L) ) 

Qsat=Q2-Qhex 
K=4 

The  pinch  point  location  is  at  the  water  entry  point  and  the 
steam  exits  the  HRSG  in  saturation. 


MDIF 
GO  TO  90 
MDIF 

last  choice,  steam  reaches  superheat  with  pinch  point  at  water 
entry  point,  entry  to  saturation  region,  or  at  steam  exit 
(approach  temperatures). 

With  assumption  pinch  point  is  at  approach  temperatures  check  to 
see  if  temperature  difference  at  other  possible  pinch  point  loca¬ 
tions  is  more  than  minimum  limit. 

T5s=T9-26.0 

CALL  CPH20(Tsat(L),T5s,Cps) 

Q3=Mr*(Cps*(T5s-Tsat(L) )+Hfg(L) ) 

Initial  guess  on  Tpp2  to  determine  Cpmix. 

Tppl =0.0 

Tpp2=Tsat ( L )+1 00 . 0 
CALL  CPM(Tpp2,T9,Mr,Cpmix) 

Tpp2=T9-Q3/( ( .95+Mr)*Cpmix) 

IF(ABS(Tpp2-Tpp1  ).GE.0.1 )  THEN 
Tppl =Tpp2 
GO  TO  30 
MDIF 
Tpp=Tpp2 

Qhex=Q3+Q1 
Tstakl =0.0 

Tstak2=T9-Qhex/ ( Cpmi x* ( Mr+0 . 95 ) ) 

IF(ABS(Tstak2-Tstak1  ).GE.0.1 )  THM 
CALL  CIW(Tstak2,T9,Mr, Cpmix) 

Tstakl  =Tstak2 
GO  TO  40 
ENDIF 

Tstack=Tstak2 

Cph=Cpmix 

CALL  CPH20(T5s,T6,Cpsc) 

Qhs=Mr*Cpec* ( T6-T5s ) 

K=5 

The  pinch  point  location  is  at  the  approach  point  temperature 
location  and  the  steam  exits  the  HRSG  as  a  superheated  vapor. 

Check  to  see  if  Mode  6  or  7  are  more  applicable. 

IF((Tstack-T1 ).LT. 26.0. AND. (Tstack-TI ).LT. (Tpp-Tsat(L) ) )  GO  TO  70 
IF( (Tpp-Tsat(L) ).LT.26.0)  THM 
Tpp=Tsat ( L ) +26 . 0 
Tstack=Tstak3 
Cph=Cpmix2 
Qhex=Ch2* ( T9-Tstack ) 

T5s2=Tsat(L)+100.0 


CALL  CPH20(Tsat(L),T5s2,Cps) 

T5s  2=  ( Qhex-Q2 )  /  (Mr*Cps  )+Tsat  ( L) 

IF(ABS(T5s2-T5s1 ).GE.0.1 )  THEN 
T5s1 =T5s2 
GO  TO  60 
ENDIF 
T5s=T5s2 

CALL  CFH20(T5s,T6,Cpsc) 

Qhs=Mr*Cpsc*(T6-T5sj 

IF( (Tstack-TI ).LT.26.0)  THEN 
GO  TO  70 
ENDIF 

K=6 

The  pinch  point  location  is  at  the  beginning  of  the  sat¬ 
uration  region  and  the  steam  exits  as  a  superheated  vapor 

GO  TO  90 
ENDIF 

GO  TO  90 

Tstack=T1+26.0 

Check  to  see  if  Mode  6  is  more  applicable. 

Tpp1=0.0 

Tpp2=Tsat ( L) +30 . 0 

CALL  CFM(Tstackf Tpp2,Mr,Cpmix) 

Tpp2=Q1  / ( (Mr+0.95 )*Cpmix)+Tstaok 
IF(ABS(Tpp2-Tpp1  J.GE.0.1 )  THEN 
Tppl =Tpp2 
GO  TO  75 
ENDIF 

IF( (Tpp2-Tsat(L) ) .LT.26.0)  THEN 
GO  TO  50 
ENDIF 

CALL  CPM(Tstack,T9,Mr,Cph) 

Qhex= (Mr+O. 95 )*Cph* (T9-Tstack ) 

T5s1=0.0 

T5s2=T9-50.0 

CALL  CFH20(Tsat(L),T5s2,Cps) 

T5s2= ( Qhex-Q2 ) / ( Mr*Cps ) +Tsat ( L ) 

IF(ABS(T5s2-T5s1 J.GE.0.1 )  THEN 
T5s1 =T5s2 
GO  TO  80 
ENDIF 
T5s=T5s2 

CALL  CFH20(T5s,T6,Cpsc) 

Qhs=Mr *Cpec* ( T6-T5s ) 


The  pinch  point  location  is  at  the  water  entry  point  and  the 
steam  exits  the  HRSG  as  a  superheated  vapor. 


RETURN 

END 


SUBROUTINE  HRSG2(T1 ,T6,T9,Mr,L,T5s,Tstack,Wes,K1  ,Qhex,Bs1  ,Bs2) 

This  subroutine  calculates  the  performance  of  a  HRSG  combined 
with  a  GT  cycle  to  produce  process  steam.  In  this  case  the  hot 
exhaust  gas  is  air  only  as  there  is  no  steam  injection.  This 
subroutine  does  the  same  calculations  as  the  HRSG1  subroutine. 

INTEGER  K1  ,L 

REAL  T1 ,T6,T9,Tpp,Tpp1 ,Tpp2,Tstack,Tstak1 ,Tstak2,Tstak3,T5s, 
+T5s1  ,T5s2,Wes(3),Cpa,Cpa1  ,Cpa2,Cps,Cps1 ,Ch1 ,Ch2,Q1  ,Q2,Q3,Qhex, 
+Qhex1  ,Tsat(5),Hfg(5),Cpw(5),Sfg(5),Mr,Bs1 ,Bs2,X 

DATA  Tstakl/0.0/,Tppl/0.0/,T5sl/0.0/,Tsat(1 )/41 7-3/, 
+Tsat(2)/444.1/,Tsat(3)/46l.7/,Tsat(4)/475-2/,Tsat(5)/486.2/, 
4flfg(l)/2132.4/,Hfg(2)/2046.2/,Hfg(3)/l984.0/,Hfg(4)/1932.7/, 
4Hfg(5)/l887.8/,Cpw(l)/4.205/,CpwT2)/4.215/,Cpw(3)/4.225/, 
-fOpw(4)/4.230/,Cpw(5)/4.245/,Sfg(1  ) /5  •  1 100/,  Sfg(2)/4. 6067/, 

4Sf g( 3 ) /4 . 2964/ , Sf g( 4 ) /4 • 0670/ , Sf g( 5 ) /3 • 8825/ 

CALL  CPH20(Tsat(L) ,T6,Cps) 

Bsl  =Mr*(Cpw(L)*(Tsat(L)-T1  *(1  .O+L0G(Tsat(L)/T1 ) ) )+ 

+Hfg(L)-T1  *Sfg(L)+Cps*(T6-Tsat  (L)-T1  *L0G(T6/Tsat  (L) ) ) ) 

Tpp=T 1+26.0 

Q1  =Mr*Cpw(L)*(Tsat(L)-T1 ) 

CALL  CPAIR1  (Tpp,T9,Cpa) 

Cpal =Cpa 

Qhexl  =0.95*Cpa1 *(T9-Tpp) 

CALL  CPH20(Tsat(L),T6,Cps1  ) 

IF ( Qhexl .LT.Q1 )  THEN 
Tstack=Tpp 

T5s=Qhex1 / (Mr*Cpw(L) )+T1 

Wes ( 1 ) =Mr*Cpw ( L ) * ( Tsat ( L ) -T5s-T1 *L0G ( Tsat ( L ) /T5s ) ) 
Wes(2)=Mr*(1 .0-T1  /Tsat(L) )*Hfg(L) 

Wes(3)=Mr*Cps1 *(T6-Tsat(L)-T1 *L0G(T6/Tsat(L) ) ) 
Bs2=Mr*Cpw(L)*(T5s-T1*(1  .Q+L0G(T5s/T1  ) )) 

K1  =1 

IF( (T9-T5s) .LT.26.0)  THEN 
T5s=T9-26.0 

Qhex=Mr*Cpw(L)*(T5s-T1 ) 

IF(Qhex.GT.Q1  )  THEN 

T5s=Qhex1 /(Mr*Cpw(L) )+T1 
GO  TO  20 
ENDIF 

Tstak2=T 1+50.0 

CALL  CPAIR1 (Tstak2,T9,Cpa) 

Tstak2=T9-Qhex/ ( 0 . 95*Cpa ) 


IF(ABS(Tstak2-Tstak1 ).GE.0.1 )  THEN 
Tstakl =Tstak2 
GO  TO  10 
ENDIF 

Ts“tack=Ts'tak2 

Wes(1 )=Mr*Cpw(L)*(Tsat(L)-T5s-T1*LOG(Tsat(L)/T5s) ) 
Wes(2)=Mr*(1.0-T1  /Tsat(L))*Hfg(L) 

Wes(3)=Mr*Cps1 *(T6-Tsat (L)-T1 *L0G(T6/Tsat (L) ) ) 
Bs2=Mr*Cpw(L)*(T5s-T1  *(1 .0+LOG(T5s/T1 ) ) ) 

K1  =2 

GO  TO  130 
ENDIF 
Qhex=Qhex1 
GO  TO  130 
ENDIF 

Tpp=Tsat ( L ) +26 . 0 

fPo  +  olfl  — O  O 

Tstak2=Tpp-Q1  /(0.95*Cpa) 

IF(ABS(Tstak2-Tstak1  ).GE.0.1 )  THEN 
CALL  CPAIR1 (Tstak2,Tpp,Cpa) 

Tstakl  =Tstak2 
GO  TO  40 
ENDIF 

Tstak3=Tstak2 

CALL  CPAIR1  (Tstak2,T9,Cpa2) 

Chi  =Cpa1  *0. 95 
Ch2=Cpa2*0 . 95 
Qhex=Ch2* (T9-Tstak2 ) 

Q2=Mr*(Hfg(L)+Cpw(L)*(Tsat(L)-T1 )) 

IF(Qhex.LT.Q2.QR.Qhex1  .LT.Q2)  THEN 
T5s=Tsat(L) 

Wes ( 2 ) = ( 1 • 0-T1 /Teat (L) ) * ( Q2-Qhex ) 

Wes(3)=Mr*Cps1  *(T6-T5s-T1  *L0G(T6/T5s) ) 

X= ( Qhex-QI ) / ( Mr*Hf g( L ) ) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1 .0+LOG(Tsat(L)/T1 )))- 
+T1  *X*Sfg(L) )+Qhex-Q1 
K1  =3 

IF((Tstack-T1 ).LT.26.0)  THEN 
Tstack=T1 +26 . 0 
Tpp2=Tsat ( L ) +30 . 0 
CALL  CPAIR1 (Tstack,Tpp2,Cpa) 

Tpp2=Q1 /(0.95*Cpa)+Tstack 
IF(ABS(Tpp2-Tpp1  ).GE.0.1 )  THEN 
Tppl =Tpp2 
GO  TO  50 
ENDIF 

IF( (Tpp2-Tsat(L) ) .LF.26.0)  GO  TO  30 
Qhex=Ch1 *(T9-Tstack) 

Wes(2)=(1 .0-T1 /Tsat(L))*(Q2-Qhex) 

Wes(3)=Mr*Cps1 *(T6-T5s-T1  *LOG(T6/T5s) ) 

X=( Qhex-QI  )/(Mr*Hfg(L)) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1 .0+L0G(Tsat(L)/T1  )))- 
+T1  *X*Sfg(L) )+Qhex-Q1 


GO  TO  130 
MDIF 

T5s=T9-26.0 

CALL  CPH20(Tsat(L),T5s,Cps) 
Q3=Mr*(Cps*(T5s-Tsat(L))+Hfg(L)) 

Tpp1=0.0 

Tpp2=Tsat ( L ) +1 00 . 0 
CALL  CPAIR1 (Tpp2,T9»Cpa) 

Tpp2=T9-Q3/(0.95*Cpa) 

IF(ABS(Tpp2-Tpp1 J.SE.0.1 )  THEN 
Tppl =Tpp2 
GO  TO  60 
ENDIF 
Tpp=Tpp2 
Qhex=Q3+Q1 
Tstak1=0.0 

Tstak2=T9-Qhex/ ( Cpa*0 . 95 ) 

IF(ABS(Tstak2-Tstak1 j.GE.0.1 )  THEN 
CALL  CPAIR1 (Tstak2,T9,Cpa) 

Tstakl =Tstak2 
GO  TO  70 
INDIP 

Tstack=Tstak2 

CALL  CFH20(T5s,T6,Cps) 

Wes(3)=Mr*Cps*(T6-T5s-T1 *L0G(T6/T5s) ) 

CALL  CPH20(T5s, Tsat(L) ,Cps) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1  .OfLOG(Tsat(L)/T1 )))+ 
•+Hfg(L)-T1  *Sfg(L)40ps*(T5s-Tsat(L)-T1 *L0G(T5s/Tsat(L) ) ) ) 
K1  =5 

IF( (Tstack-TI ). LT. 26.0. AND. (Tstack-TI ).LT. (Tpp-Tsat(L) ) ) 
-KX)  TO  100 

IF( ( Tpp-Tsat (L) ) . LT . 26 . 0)  THEN 
Tpp=Tsat(L)+26. 0 
Tstack=Tstak3 
Qhex=Ch2* ( T9-Tstack ) 

T5s2=Tsat (L)+1 00. 0 

CALL  CPH20(Tsat(L),T5s2,Cps) 

T5s2=(Qhex-Q2)/(Mr*Cps)+Tsat(l) 

IF(ABS(T5s2-T5s1  J.GE.0.1 )  THEN 
T5s1 =T5s2 
GO  TO  90 
MDIF 
T5s=T5s2 

CALL  CPH20(T5s,T6,Cps) 

¥ea(3)=Mr*Cps*(T6-T5s-T1  *L0G(T6/T5s) ) 

CALL  CPH20(Tsat{L),T5s.Cpe) 

Bs2=Mr* ( Cpw ( L) * ( Tsat ( L)-T1  * ( 1  .OfLOG(Tsat(L)/T1  )))+ 
4Hfg(L)-T1 *Sfg(L)4Cpe*(T5s-Tsat (L)-T1 *L0G(T5s/Tsat (L) ) ) ) 
K1=6 

IF( (Tstack-TI ).LT. 26.0)  GO  TO  100 
GO  TO  130 
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GO  TO  130 

1 00  Tstack=T  1+26.0 

Tpp1=0.0 

Tpp2=Tsat(L)+30.0 
110  CALL  CPAIR1  (Tstack,Tpp2,Cpa) 

Tpp2=Q1 / (0.95*Cpa)+Tstack 
IF(ABS(Tpp2-Tpp1 J.GE.0.1 )  THEN 
Tppl =Tpp2 
GO  TO  110 
ENDIF 

IF((Tpp2-Tsat(L)).LT.26.0)  GO  TO  80 
Qhex=Ch1  *(T9-Tstack) 

T5s1=0.0 

T5s2=T9-50.0 

120  CALL  CPH20(Tsat(L),T5s2,Cps) 

T5s2= ( Qhex-Q2 ) / (Mr*Cps ) +Tsat ( L ) 

IF(ABS(T5s2-T5s1 ).GE.0.1 )  THEN 
T5s1=T5s2 
GO  TO  120 
ENDIF 
T5s=T5s2 

CALL  CPH20(T5s,T6,Cps) 

Wes ( 3 ) =Mr*Cps* ( T6-T5s-T1 *L0G(T6/T5s) ) 

CALL  CPH20(Tsat(L),T5s,Cps) 

Bs2=Mr*(Cpw(L)*(Tsat(L)-T1*(1 .O+L0G(Tsat(L)/T1 ) ) )+ 
+Hfg(L)-T1  *Sfg(L)+Ops*(T5s-Tsat(I)-T1  *L0G(T5s/Tsat  (L) ) ) ) 
K1  =7 

130  RETURN 
END 


SUBROUTINE  CPAIR1 (T1 ,T2,Cpa) 

This  subroutine  finds  the  specific  heat  at  constant  pressure 
of  air  as  a  function  of  average  temperature  of  the  process. 

It  is  based  on  a  relation  by  Irvine  and  Li ley,  STEAM  AND  GAS 
TABLES  WITH  COMPUTER  EQUATIONS,  Academic  Press,  Inc.,  Orlando, 
EL,  1984. 

INTEGER  I 

REAL  T1,T2,Tavg,A(5),Cpa 
C 

DATA  A(1 )/0.103409E1 /,A(2)/-0.2848870E-3/,A(3)/0.781 6818E-6/, 
+A ( 4 ) /-0 • 4970786E-9/ , A ( 5 ) /0 . 1 077024E-1 2/ 

C 

Cpa=0.0 

Tavg= (T1 +T2 )/2 . 0 
DO  10  1=1,5 

Cpa=Cpa+A ( I ) * ( Tavg** ( 1-1 )) 

10  CONTINUE 
RETURN 
END 


oooo  ooo  o  o  ooooo  ooo  -*■  o  o  ooooooo 


SUBROUTINE  CHi20(T1 ,T2,Cps) 

This  subroutine  finds  the  specific  heat  at  constant  pressure 
of  steam  (ideal  gas)  as  a  function  of  average  temperature  of 
the  process.  It  is  based  on  a  relation  by  Reynolds,  THERMO¬ 
DYNAMIC  PROPERTIES  IN  SI,  Stanford  University,  Stanford,  CA, 
1979- 

INTEGER  I 

REAL  T1 ,T2,Tavg,G(6),Rs,Cvs,Cps 

DATA  G(1  )/4.6E+04/,G(2)/l  .01 1 249*>+03/,G(3)/8.3893E-01  /, 

-KJ  ( 4 ) /-2 . 1 9989B-04 / , G ( 5 ) /2 . 466 1 9E-07 / , G ( 6 ) /-9 . 7047E- 1 1 / , 
+Rs/.461 52/ 

Cvs=0.0 

Tavg=(T1+T2)/2.0 
DO  10  1=1,6 

Cvs=Cvs+G ( I ) * ( Tavg** ( 1-2 ) ) 

CONTINUE 

Cps=(Cvs/l000. )+Rs 

RETURN 

END 


SUBROUTINE  CiW(T1  ,T2,Mr,Cpmix) 

This  subroutine  finds  the  specific  heat  at  constant  pressure 
of  the  steam/air  mixture  of  the  STIG  cycle  using  the  two 
subroutines  listed  above. 

REAL  T1 ,T2,Mr,Cpa,Cps,Cpmix 

CALL  CPAIR1 (T1 ,T2,Cpa) 

CALL  CPH20(T1 ,T2,Cps) 

Cpmix=  (Mr*Cps-K) .  95*Cpa )  /  ( 0 . 95+Mr ) 

RETURN 

END 

- REAL  FUNCTIONS - 

REAL  FUNCTION  Tisen(TYPE,T1 ,Pr,Ka) 

This  function  finds  the  end  state  temperature  for  an  isentropic 
compression  or  expansion  process. 

CHARACTERS  TYPE 
REAL  T1,Pr,Ka 


IF(TYPE.EQ.'COMP’)  THEN 


Tisen=T1*(Pr**((Ka-1 .0)/Ka)) 
TiTSP 

Tisen=T1 /(Pr**( (Ka-1 .0)/Ka) ) 
ENDIP 
RETURN 
END 


REAL  FUNCTION  Tact ( TYPE, T1 , Ts2 , ETA ) 

This  function  finds  the  actual  end  state  temperature  for  a 
compression  or  expansion  process  using  the  isentropic 
process  end  temperature  and  the  component's  isentropic 
efficiency. 

CHARACTERS  TYPE 
REAL  T1 ,Ts2,ETA 

IF(TYPE.EQ. 'COMP' )  THEN 
Tact=T1 + ( Ts2-T1 ) /ETA 
USE 

Tact=T1-ETA*(T1-Ts2) 

ENDIP 

RETURN 
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PROGRAM  CONDEN 

This  program  modifies  data  files  to  insure  no  steam  condenses 
inside  the  components  of  the  STIG  cycle.  It  checks  the  output 
files  of  max  First  Law  efficiency  and  corresponding  steam/air 
mass  ratio  versus  turbine  inlet  temperature.  If  the  stack  temp 
at  max  efficiency  is  below  the  dewpoint  temp  then  the  mass  ratio 
is  decreased  until  the  stack  temp  is  above  dewpoint.  This  program 
must  be  linked  with  the  subroutines  of  PROGRAM  GT. 

CHARACTER*!  CYCLE 
CHARACTERS  3  FNAME 
INTEGER  I, J,K,K1  ,L,NN,SIZE 

REAL  MAX(1334,5),T(9),T6,Tstack,Tstak1 ,T5s,T5s1 ,Tdp,T9a,P(9),PL, 
+ETA,EPS(6),Wc(2),Wt(2),Wp,Wnet,Cph,Cpsc,Cpmix,Qha(2) ,Qhs,Qhex, 
+Qhexa,Be,We(3),Wt1 ,Bia,Biw,Be1 ,Wes(3) ,Bs1 ,Bs2,Mr 

CYCLERS' 

WRITE(*,10) 

FORMAT ( '0' ,  'What  is  the  file  to  be  reviewed  (FOR  .DAT;  )?’) 
READ(*,20)  FNAME 
FORMAT (A1 3) 

OPPN  (UNIT=1 0 , NAME=FHAME , STATUS^ ' UNKNOWN • ) 


NN=1 

READ  (10,*,END=40)  (MAX(NN,J),  J=1,5),L 

NN=NN+1 

GO  TO  30 

CONTINUE 

CLOSE (UNIT=1 0,STATUS='SAVE' ) 

SIZE=NN-1 
DO  80  1=1, SIZE 
T6=MAX(I,1 ) 

Mr=MAX(I,2) 

Tstack=MAX( 1,5) 

CALL  DEWPT(Mr,Tdp) 

IF(Tstack.LI.Tdp)  THEN 


DO  50  J=1 ,500 
Mr =Mr -0.001 
CALL  DEWPT(Mr,Tdp) 

CALL  CYCLE2(T6,L,Mr,CYCLE,T,P,PL,ETA,EPS,Wc,Wt,Wp, 
+Wnet,Cph,Cpsc,Cpmix,Qha,Qhs,Qhex,Tstack,T5s,K,Be,We,Wt1  ,Bia, 
+Biw,Be1 ,Wes,K1 ,T5s1 ,Tstak1 ,Qhexa,T9a,Bs1  ,Bs2) 
IF(Tstack.GT.Tdp)  GO  TO  60 
CONTINUE 

MAX(I,2)=Mr 

MAX(I,3)=ETA 


MAX(I,4)=Wnet 

MAX(I,5)=Tstack 

ENDIF 

WRITE(19,70)  (MAX(I,J),  J=1 ,5) 

FORMAT ('  ,,F6.1,2X,F7.3,2X,P6.4,kX,F8.3,2X,F9.5) 
CONTINUE 

STOP 

END 


■SUBROUTINE 


SUBROUTINE  DEWPT(Mr,Tdp) 

This  subroutine  determines  steam  saturation  temperature  as  a 
function  of  saturation  pressure.  It  is  based  on  a  relation 
by  Irvine  and  Liley,  STEAM  AND  GAS  TABLES  WITH  COMPUTER 
EQUATIONS,  Academic  Press,  Inc.,  Orlando,  PL,  1984. 

REAL  Mr,Psat,Ys,Tdp 

DATA  MWa/28 . 97 / , MWs / 1 8 . 1 01 5 / , A/42 . 6776 / , B/-3892 . 70 /, C /-9 . 48654 / 

Ys= ( Mr*MWa ) / ( 0 . 95*MWs+Mr*MWa ) 

Psat=Ys*1 01 .325 

Tdp=A+B/( ALOG ( Psat/1 000 . 0)4C ) 


RETURN 
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